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ABSTRACT 


hodern  U.S.  military  radars  and  communicat  ion  equipment 
performance  will  be  strongly  influenced  by  the  environment 
it  will  be  operating  in.  One  of  tie  most  important  atmos¬ 
pheric  affects  is  ducting  of  electro  mag  ne  tic  energy  by  re¬ 
fractive  layers  in  the  atmosphere.  To  assess  the  affects  of 
ducting  on  electromagnetic  emissions  around  Dhanran,  Saudi 
Arabia,  a  geometric  optics  model  of  wave  propagation  devel¬ 
oped  by  Raymond  P.  Was!cy  was  modified  and  utilized.  This 
thesis  also  attempted  to  show  the  correlation  of  wind  direc¬ 
ti  cn  and  wind  speed  to  the  establishment,  location  and 
heignts  of  lucts.  Finally  this  thesis  attempted  no  deter¬ 
mine  if  there  was  any  correlation  between  the  occurrence  of 
land-sea  breezes  and  ducting. 
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I. 


IHTR  CCDCTIO N 


A.  BACKGROUND 

The  United  States  (U.S.)  military  anticipates  fighting 
outnumber ea  in  ar.y  future  conflict.  Since  Vietnam,  modern¬ 
ization  of  our  forces  has  dragged,  thereby  permitting  a  de¬ 
cay  in  the  readiness  cf  our  forces.  At  the  same  tame,  the 
Soviet  Union  has  increased  defense  expenditures.  The  Sovi¬ 
ets  remain  committed  to  their  goal  of  world  socialism  and 
Soviet  Policy  proceeds  on  the  lasis  of  military  power. 


— 

Country  1973  1  979  1  93  1 

United  States  78,472  114,503  171,000 
Soviet  Union  92,000  241,000  2o7,U00 

Figure  1;  ’Defense 
Expenditures  (thousand  $) 

[Ref.  1,  2,  S  3]  « 

Major  Army  unit  commanders  must  be  able  to  detect,  track 
and  destroy  enemy  targets  deep  in  the  enemy  territory  before 
they  reach  the  f crwar d-line-of -cwn-t r cop s  (FLOT) .  Our  com¬ 
manders  must  have  the  time  to  adjust  their  own  forces  to  be 
able  to  meet  the  enemy  on  a  mere  even  combat  ratio  when  the 
enemy  reaches  the  ELOT.  Thru  an  integrated  collection  ef¬ 
fort  and  the  use  of  sophisticated  electronic  equipment,  it 
is  the  jot  of  the  intelligence  community  to  find  and  locate 
the  enemy.  Tc  accomplish  the  above,  effective  use  of  the 
electromagnetic  spectrum  is  essential. 

The  performance  of  any  electromagnetic  system  varies  with 
the  geographical  region,  tepegraphy,  systems  employment 
configuration  and  operating  frequency.  These  performance 


however  m 
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character istics  are  recognized  by  most  planners, 
addition  to  these  parameters  one  must  also  consider  meteoro¬ 
logical  elements.  Most  testing  of  electronic  equipment  is 
performed  ander  normal  conditions.  This  is  a  tactical 
shortcoming  because  extreme  environmental  conditions  may 
have  an  adverse  affect  on  a  systems  performance. 

B.  OBJECTIVE 

The  objective  ox  this  thesis  is  to  determine  the  synoptic 
meteoro  logical  conditions  that  will  severely  affect  electro¬ 
magnetic  (EM)  propagation  on  the  Arabian  Peninsula.  Atmos¬ 
pheric  refracti vity,  surface  and  elevated  ducts,  and  land- 
sea  breeze  will  be  addressed. 

C.  THE  DESERT  ENVIRONMENT 

Tne  desert  environment  was  selected  as  a  type  of  environ- 
mental/topc graphical  region  for  examining  the  affects  of  ex¬ 
tremes.  Radiosonde  data  recorded  in  Dhahran,  Saudi  Arabia 
from  1978  to  1  980  w.ill  be  examined.  The  importance  of  this 
data  is  for  several  reasons.  Eirst,  the  desert  areas  of  tne 
world  comprise  approximately  191?  of  the  earth’s  land  mass. 
This  is  a  significant  portion  of  the  total  surface  area 
available  for  ground  combat.  Saudi  Arabia  encompasses  an 
area  of  330,000  sguare  miles,  much  of  which  is  desert  [Ref. 
4].  Second,  several  of  the  desert  regions  of  tne  world  have 
a  political  and  military  significance  because  of  their  stra¬ 
tegic  location  and  valuable  mineral  deposits.  The  Arabian 
oil  fields  account  for  about  cne-haif  of  tne  known  reserves 
of  the  non-Co  mm un ist  world.  They  also  supply  atout  one- 
fifth  of  the  world’s  total  oil  production  [Ref.  5].  Saudi 
Arabia  was  the  largest  producer  of  crude  petroleum  in 
the  Middle  East  and  the  third  largest  in  the  world  in 


Ccun  try 

Soviet  Union 
United  States 
Saudi  Arabia 
Iran 

Venezula 

Kuwait 

Ira.-, 


1974 

3 , 373 ,0  50 
3,  202,565 
2,596, 543 
2,  197,  700 
1,086,333 
394, 781 
720, 729 


197  9 

4,307,  100 
3,111,  625 
3,479, 339 
1,121, 346 
860, 072 
602, 000 
1,252, 000 


1  980 

4,432, 050 
94  7,  905 
3, 530, 000 
550, 000 
793,397 
912, 610 
361,000 


Figure  2:  'Crude 
Oil  Production  Comparisons 
(thousand  42-gal  barrels) 


[  Ref. 


7,8, 


9] 


- 1 


1980  [Ref.  10].  Thirdly,  the  desert  has  the  most  extreme 
cases  of  abnormal  atmospheric  ref rac tivit y .  Fourthly,  a 
worldwide  analysis  of  upper  atmospheric  radiosonde  data  was 
performed  by  Crtenturger  (GTE,  Sylvania  1973)  .  Results  in¬ 
dicate  that  trade  wind  regions  are  the  areas  of  significant 
ducting.  One  area  of  prevalent  ducting  was  found  in  the  up¬ 
per  Indian  Ocean  and  Persian  Gulf.  The  probability  of  duct¬ 
ing  in  this  area  is  6  0%  [Ref.  11].  Finally,  the  dramatic 
increase  of  Soviet  military  power  in  Asia,  the  Pacific  Ocean 
and  the  Indian  Ocean  is  the  most  significant  military  devel¬ 
opment  hn  recent  years.  The  Soviets  have  replaced  and  up¬ 
graded  eguipment  on  their  Southern  Theater.  Additionally , 
their  105,000  troops  in  Afghanistan  have  altered  the  balance 
of  power  in  the  Indian  Ocean  ar.d  Persian  Gulf  [Ref.  12]. 
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II.  GENERAL  discussion,  inters  EIATIQN5HIPS  and 

EQUATIONS 


A.  REFRACT! V  III 


I he  transmission  of  electromagnetic  (HM)  signals  tnrough 
a  medium  is  affected  by  the  absorption  and  re-emission  of  Ed 
energy  by  the  atomic  and  molecular  elements  of  that  medium. 
The  dielectric  constant  (£)  best  describes  the  ii.teraction 
of  the  electric  field  with  the  medium.  As  the  Ed  wave  in¬ 
teracts  with  the  new  medium,  its  speed  changes  and  is  deter¬ 


mined  by: 


C 


( £  Q  .vj  i  j 


where  c  is  the  velocity  in  the  vacuum  and  v  is  the  velocity 
in  the  medium.  Rather  than  deal  with  a  velocity,  physicists 
defined  a  new  parameter  "Index  cf  Refraction"  ( n)  wnere 

C  t - - 

n  —  -\J(c  ^  E  g’ m  2.  ,) 


A  measure  of  the  amount  of  refraction  experienced  Ly  a 
ray  as  it  passes  through  a  surface  which  separates  two  media 
of  different  densities  is  designated  "n".  It  is  tre  ratio 
of  the  wavelength  cr  velocity  cf  an  EM  wave  in  a  vacuum  to 
that  in  the  new  medium.  When  the  EM  wave  passes  from  one 
medium  to  another  non-absorbing  medium,  the  angle  of  inci¬ 
dence  e,  and  the  angle  of  refraction  02  (See  Figure  3)  are  re¬ 
lated  by  the  principles  of  Snell’s  Law: 

9.‘  -  '2±.  (tgM  3) 

vn  6,  n  , 

I  be  Index  cf  Refraction  of  a  vacuum  is  one;  of  air  is  ap¬ 
proximately  1.0  00326;  and,  for  water  approximately  1.33. 
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n2>a1 

V°1 


I  Figure  3 : 


'Refraction  of  Incident  Rays 
'Ref.  13]  • 


The  relationship  of  the  index  of  refraction  to  the  atmos¬ 
pheric  pressure  (?)  ,  water  vapor  pressure  (e)  ,  and  teapera- 
ture  (T)  is  given  ly  the  following  equation: 


Cn-0  -  (77  <« 


5".  1  3  *■  3.  7  s"  *  l  C  3  — -  )  x  j  0 

i  -rz 


(ig  m  H) 


where  atmospheric  pressure  and  the  water  vapor  pressure  are 
ir.  mrllihars  and  temperature  is  in  degrees  Kelvin.  The  re¬ 
sulting  values  of  the  index  of  refraction  are  aw  h  war  d . 
Therefore  a  new  parameter  "N",  which  is  the  ref rac ti vi t y,  is 
defined  for  convenience  as  follows: 


M  -  (  r\  -  >  )  *■  I  C,  . 


C  £g  u  T  ) 


*.  :i  us , 


N  -  77.4  zr 


e  ,  -  ,  .>  e. 

-  S  L  —  V  .-5  7o>  *  /C  — 

r  t 


( 6  ) 


f  R  e  f .  la]. 
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B.  MODIFIED  REFRACTIVITY 

The  use  of  refractivity  to  depict  refraction  is  diffi¬ 
cult.  Ducting  conditions  are  identified  by  determining  its 
char.je  with  height.  In  a  normal  situation,  the  r  efr  act  iv  i  t  y 
decreases  with  height.  The  refractivity  nay  be  aathexati - 
cally  modified  so  that  when  its  gradient  (dN/dZ)  is  applied 
to  EM  propagation  over  a  hypothetical  flat  eartn  it  is  es¬ 
sentially  equivalent  to  the  propagation  over  the  real  curved 
earth  with  the  actual  refractivity. 

The  Modified  Refractive  Index,  M,  adds  157  N-ur.its  rer 
kilometer  to  all  N-values.  It  is  defined  as: 

|V(  -  N  r  *  ^  7  c.  •  t  s  m  '<  ) 

where  N  is  the  value  of  refractivity  at  any  height  Z,  and  Z 
is  un  Kilometers.  M  will  increase  with  height  in  the  stan¬ 
dard  atmosphere.  Also,  when  the  M-gradient  (dM/dZJ  is  zero 
the  ray  curvature  is  egual  to  the  real  earth's  curvature.  M 
versus  height  profiles  are  mainly  used  to  obtain  ducting  in¬ 
formation  (See  Chapter  II,  Section  D)  .  [Ref.  15] 

C.  DN/DZ  AND  D  M/D  Z 

An  important  description  of  refractivity  is  not  its  value 
but  its  gradient.  Refractivity  is  a  multi-variable  parame¬ 
ter  and  thus  an  expression  for  its  gradient  is: 

^Kl  \L  ‘I-  5;  ) 

£  ^  ■}  P  A1  i  ^  1  ^  ^  ^  cV  * 

For  an  average  standard  condition,  dN/dZ  is  approximately 
-40  per  Km.  In  ray  tracing  problems,  the  gradient  dM/dZ  can 
be  use-  i  to  ortain  a  ray  path  curvature  fnat  is  relative  to 
the  curvatULe  cf  the  earth. 


Clas  si f icat ion 


uK/dZ  (/Km)  dt</dZ  (/Km)  Range 


Normal 

Super  Refraction 
Trap  pi ng 


>0 

to  - 7 9 
tc  -157 
<-  157 


>157 

79  to  157 
0  to  79 
<0 


Reduced 
Nor  mal 
Increased 
m creased 


Figure  4:  ' N  £  M  Gradients 

[Ref.  16]  ’ 


_ j 


SuDfoiraction  ^  0 


Figure  5:  'N  Gradient  Pictorial 
[Ref.  17]  ’ 


Standard  (normal)  propagation  results  in  a  ray  curvature 
due  to  refraction  which  has  a  value  approximately  one-fourth 
that  of  the  earth's  curvature.  This  is  equivalent  to  the 
straight  line  propagation  over  a  hypothetical  earth  whose 
radius  is  four-thirds  the  radius  of  actual  earth.  Subre¬ 
fraction  produces  a  less  than  normal  downward  bending  or 


ever,  an  upward  tending  of  radio  waves  as  they  travel  through 
the  atmosphere.  Thus,  radar  and  radio  coverage  is  le- 
creased.  Super  refraction  produces  a  greater  than  normal 
downward  tending  of  radio  waves  as  they  travel  through  the 
atmosphere.  This  results  in  extended  radio  horizons  and  in¬ 
creased  radar  coverage.  Strong  super  refraction  can  produce 
skip  effects  in  elevated  layers  or  the  troposphere.  Skip 
effects  occasionally  rake  it  possible  to  detect  targets  at 
distances  greater  than  the  normal  horizon  while  closer  tar¬ 
gets  remain  undetected. 

E.  DUCTING 

A  duct  is  a  shallow,  almost  horizontal  layer  in  the  at¬ 
mosphere  where  EM  energy  is  trapped.  The  trapping  layer  is 
where  dN/dZ  <  -157  Km  and  thus  the  ray  will  be  bent  toward 
the  earth.  The  trapping  layer  is  the  top  of  the  duct.  En¬ 
ergy  transmitted  within  the  duct  will  be  partially  confined 
and  channeled  between  the  top  and  bottom,  of  the  duct.  Duct¬ 
ing  occurs  in  several  ways  and  can  best  be  categorized  by 
the  altitudes  at  which  they  are  found: 

Surface  Ducts  -average  height  less  than  1,500  feet; 

Elevated  D  uct  s-a  vera  <je  height  5,OC0  to  10,000  feet; 

Evaporation  E  ucts-approxi  irately  100  feet. 

Evaporation  ducts  cnly  occur  over  water. 

These  duct  types  are  defined  by  Figures  6  thru  8.  The 
top  of  a  duct  corresponds  to  a  height  above  the  surface 
where  the  M  value  is  a  minimum.  The  duct  base  corresponds 
to  the  height  at  which  a  vertical  line  drawn  downward  from 
the  point  of  minimum  value  first  intersects  a  point  of 
egual  A  units  or  the  surface.  Eadar  targets  may  be  detected 
at  long  ranges  in  both  target  and  radar  are  in  the  duct. 
The  area  just  above  a  duct  may  have  reduced  radar  coverage. 
An  aircraft  or  missile  flying  just  above  the  duct  mignt  not 
be  detected  until  very  close  to  the  radar,  if  at  all. 


Figure  7:  ’Surface 
Eased  Duct  Caused  fcy  Surface  Layer 


elevated  layer 


Figure  8:  ’Surface 

uct  Caused  ty  ar.  Elevated  Layer 
£  Pef .  20  ]  ’ 
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A  luct  Mill  trap  EM  energy  lor  only  a  selected  frequency 
ran  ye .  Minimus:  trapping  frequencies  lave  been  established 
because  the  limit  is  cn  the  lav*  side.  The  minimum  frequency 
that  will  he  trapped  is  given  fcy: 

f  _  .  .  ..  M 

r  m  - 


o.  G  0  33  -  10  J 


where  d  is  the  thickness  of  the  duct  in  meters.  [Ref.  22] 


E.  LAND-SEA  BREEZE 

The  land-sea  breeze  is  the  complete  cycle  of  the  day- 
night  local  winds  cccuring  on  sea  coasts  due  to  the  differ¬ 
ences  in  surface  temperature  of  the  land  and  sea.  The  land 
breeze  component  of  the  system  blows  from  land  to  sea  and 
the  sea  breeze  blows  from  sea  tc  land. 

The  basic  principle  of  the  sea  breeze  is  that  during  the 
day  the  land  and  the  air  over  the  land  gets  heated  consider¬ 
ably,  while  the  air  over  the  sea  changes  slightly.  The 
warm,  light  air  over  the  land  then  rises  and  is  replaced  by 
the  cooler  air  from  the  sea.  The  day  time  sea  breeze  sur¬ 
passes  in  intensity  the  night  time  land  breeze.  The  direc¬ 
tion  of  the  sea  breeze  does  not  remain  constant  during  the 
course  of  tne  day.  The  gradual  change  in  the  direction  of 
the  sea  breeze  appears  because  of  the  affect  of  the  Coriolis 
Force. 

At  night  the  wind  direction  reverses,  because  the  air 
over  the  land  becomes  cooler  than  the  sea  air.  Now  it  is 
the  sea  air  which  rises  and  the  cooler  land  air  that  moves 
out  from  the  land  as  a  land  breeze.  [Ref.  23] 


1800 
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Figure  10:  'Hourly 
Variation  of  Relative  Wine  Velocity 
at  Hcek  van  Holland  31  July  1533 
[Ref.  24  ]  • 


Ill*  SAODI  ASAEIAN  H EATHER 

A.  THE  CCONTRY 

Saudi  Arabia  occupies  approximately  four-fifths  of  tne 
Arabia;.  Peninsula  (see  Figure  11).  It  is  aoout  the  sire  of 
the  United  States  east  of  the  Mississippi.  The  Arabian  Pen¬ 
insula  is  a  plateau  which  slopes  slightly  toward  the  east. 
It  contains  both  the  world's  largest  sand  lesert,  the  Rut 
al-Khaii,  and  maybe  the  world's  largest  oasis,  al-Hasa.  In 
addition  to  the  Put  al-Khali,  cr  "Empty  Quarter",  the  otner 
two  sand  areas  are  the  Great  Nufud  Desert  and  tne  Dahana  De¬ 
sert.  Outside  these  deserts  the  surface  is  gravel,  or  in 
the  case  of  the  west-central  area  the  surface  consists  of 
crumbled  beds  of  lava.  [Ref.  25  j 

E.  CLIMATOLOGY 

Saudi  Arabia  has  a  desert  climate  characterized  by  ex¬ 
treme  heat  during  the  day  with  an  abrupt  drop  in  temperature 
at  night  and,  a  small  but  erratic  rainfall.  Along  the 
coastal  regions  of  the  Red  Sea  and  the  Persian  Gulf  the  de¬ 
sert  temperature  is  moderated  by  the  closeness  of  those  bod¬ 
ies  of  water.  Temperatures  seldom  go  ever  loo0’?,  but  the 
relative  humidity  is  usually  over  85  percent.  This  combina¬ 
tion  produces  a  hot  mist  during  the  day  and  a  warm  fog  at 
night.  In  Uajd  and  the  deserts  a  uniform  climate  prevails. 
The  average  temperature  is  112°F.  Readings  of  up  to  130^F 
are  common.  In  the  winter  the  temperature  seldom  drops  be¬ 
low  32UF,  however,  the  almost  total  absence  of  humidity  and 
high  wind-chill  factor  make  for  a  cold  atmosphere.  In 
spring  and  autumn  the  temperatures  average  85UF.  Along  the 
western  and  eastern  coastal  strips  the  prevailing  winds  are 
from  the  northern  quadrant.  A  southerly  wind  is  accompanied 
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Figure  1 2 : 

'Mean  Temperatures  and  Precip 
Values  [kef  26  £  27]  ' 

i  t  a  t  i  o  n 

by  an  increase  in  temperature  and  humidity  and  by  a  particu¬ 
lar  kind  of  storm  known  in  the  Persian  Gulf  area  as  "kauf". 
In  late  spring  and  early  summer  a  strong  northwesterly  wind 
called  a  "shamel",  blows  particularly  severe  in  eastern  Ara¬ 
bia.  The  shamel  produces  sand  and  dust  storms. 
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In  winter,  the  Mediterranean  cyclones  (lows)  moving  west 
to  east  in  association  with  upper  troughs  and  jet  streams, 
are  considered  the  main  rain  producing  systems.  The  distri¬ 
bution  of  winter  rainfall  shows  maximum  values  in  the  penin¬ 
sula's  northern  part  and  a  gradual  decrease  in  the  lowlands 
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Figure  15; 


on  ta«i  eastern  an  d  western  sides.  In  the  sprit  g,  the  doiz- 
terrar.ean  cycicr.es  continue  to  afreet  the  north.  Jyn:. g  is 
also  characterized  by  weak  stability  in  the  lower  .itaosrr. er¬ 
ic  layers  and  large  daytime  differences  between  land  and  wa¬ 
ter  surfaces.  These  conditions  stimulate  active  local  cir¬ 
culations  between  land  and  sea,  and  between  mountains  ar.  d 
valleys  particularly  in  tne  southwest.  In  summer,  tne  ther¬ 
mal  monsccn  treugn  is  established  across  the  peninsula. 
Thus  rainfall  is  restricted  tc  the  peninsula's  soutawest 
while  the  northern  part  is  dry.  In  autumn,  middle  latitude 
disturbances  begin  to  affect  the  northern  portion  or  tne 
peninsula  but  local  circulation  is  weak  therefore  there  is 
only  minor  concentrations  of  rainfall  in  the  north.  The 
frequency  of  thunderstorms  on  the  Arabian  Peninsula  is  re¬ 
lated  to  tie  rainfall  distribution.  There  is  a  higher  rre- 
guency  distribution  of  thunderstorms  near  the  eastern  and 
western  coasts  than  in  the  interior  areas.  LRef.  30] 

C.  SYNOPTIC  CCKDITICNS 

The  weather  situation  during  the  months  of  day  and  Octo¬ 
ber  for  1  973,  1  973  ,  and  1  930  in  Dnahran,  Saudi  Arabia  will 
be  overviewed  in  this  section.  Information  was  taken  from 
weather  maps  provided  by  the  Navai  Postgraduate  School  de- 
teorologicai  Department.  Data  from  radiosonde  launcaes  at 
Station  40427,  Eahrain  duharrag,  were  utilised.  3ahrair.  du¬ 
harrag  is  located  approximately  25  miles  east  of  Dhahran . 
It  snouid  also  be  noted  that  ail  watches  in  Saudi  Arabia  are 
reset  it  sundown.  The  information  posted  on  the  weather 
maps  utilized  the  Zulu  Time  Zone  System.  Saudi  Arabia  fails 
into  tne  Delta  time  zone,  thus  00  90  Z  is  04000  or  4  o'clock 
in  tne  iiornin  ; . 

Figures  15  and  17  show  the  pressure,  temperature  and  lew 
point  averages  for  1973,  1979,  and  1930.  dean  daily  temper¬ 
atures  differed  sligntly  from  history  (see  Figure  12).  day 
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had  a  mean  of  84.7'*?  (29.3<*C)  which  was  4 *  f  cooler  than  the 
norm.  October  had  a  mean  of  85.3'F  (2S.9*C)  waicn  was  2CF 
warmer  than  normal. 

Cyclones  predominantly  (45^)  affected  the  area  durir.j  the 
month  of  Kay.  Anticyclones  were  only  in  the  area  6%  of  the 
time.  however,  durir.j  Kay  197  9  anticyclones  affect  was  122. 
During  October  1979  and  1930  ,  cyclones  dominated  the  area 
30$  of  the  time.  In  October  1973,  anticyclones  were  present 


33>.  of  the  time  while  cyclones  only  2%.  Also  during  October 
haze  was  present  52 T.  cf  the  time.  Apart  iron  tne  two  excep¬ 
tions  these  percentages  are  consistent  with  tr.e  climatology 
or  the  Dhahran-  Bahrain  Muharrag  area  since  Mediterranean 
cyclones  affect  the  northern  portion  cf  tre  peninsula  durm  : 
spring  and  late  autum. 

Figures  18  and  19  sho  w  that  501  of  the  time  the  wind  s 
were  out  of  the  north  and  northwest.  Winds  out  of  the 
soutnern  guadrant  were  occasionally  seen  at  0  0  302  in  both 
May  and  October  hut  were  absent  at  12002.  winds  out  of  the 
r.ortn  and  northwest  were  the  strongest  averaging  11  knots 
(see  Figure  20). 
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Ducting  phenomenon  for  1978,  1979,  and  1980  is  presented 
in  Chapter  IV.  Fcr  comparative  purposes  results  of  two  Ra¬ 
diosonde  Data  Analysis  Projects  by  GTE  Sylvania,  Inc.  are 
included  here.  These  projects  covered  the  years  1v66  to 
1969  and  1973  tc  19  7  4.  Radiosonde  data  from  Station  4  0427, 
Bahrain  Muharraa,  provided  by  the  OSAF  Environmental  Techni¬ 
cal  Applications  Center  (ETAC)  was  utilized  by  3TE  Sylvania, 
Inc.  to  obtain  their  results.  Both  duct  and  s uper- r ef r act - 
ing  layer  (SELF.)  gradients  (d M/dZ  <  -100  N  units/Km)  are  de¬ 
picted  on  the  fcllcwing  graphs. 

Figures  21  thru  28  reflect  the  following  res;!ts  for  t:.e 
months  of  May  and  October: 

May  October 

1.  Percent  Occurrence 

Elevated  layers  10%  25% 

2.  Minimum  Trapping  Freyuency 

for  Elevated  Ducts  400  MHz  200  MHz 

3.  Percent  Occurrence 

Surface  Layers  85? 


o5% 


1  W  — *  O  y?  w  w  i+ovj 


Meters 

Figure  27:  'Thickness  Surface  Layers 
J  [  Eef .  37]  ' 
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A.  DATA  EESCEIPTICN 

Radiosonde  data  recorded  at  120QZ  (1600  Hoars  local)  and 
000QZ  (24C0  Hoars  local)  on  each  day  of  Hay  and  October  for 
1973,  1979,  and  1980  were  utilized  for  this  study.  This 
data  was  obtained  from  the  U.S.  Air  Force  Environmental 
Technical  Applications  Center  (ETAC)  at  Scott  Air  Force 
Base,  Illinois. 

This  raw  data  had  to  be  manipulated  into  a  manageable 
form  first.  Utilizing  the  Naval  Postgraduate  School  I3H 
3033  Computer  the  data  for  each  radiosonde  launch  was  organ¬ 
ized  and  cut  off  at  the  100  millibar  pressure  level.  This 
corresponded  to  a  maximum  altitude  of  approximately  16500 
meters  (54000  feet).  Next,  for  each  observation  the  refrac- 
tivity,  modified  r ef racti vit y ,  dN/dZ,  and  dM/dZ  were  comput¬ 
ed.  Each  launch  record  was  then  examined  to  determine  if 
ducting  was  present.  For  each  launch  recording  where  duct¬ 
ing  was  probanle  a  plot  of  N,  a  plot  of  H,  a  ray  trace,  and 
a  plot  of  power  density  were  dene.  Comparing  the  modified 
refractivity  plcts  and  also  the  ray  traces  nine  (9)  distinct 
groupings  were  established  for  further  analysis. 

The  computer  program  utilized  in  this  thesis  was  devel¬ 
oped  by  Raymond  P.  Nasky.  He  developed  the  model  to  ana¬ 
lyze  the  effects  of  atmospheric  refraction  on  the  field 
strength  of  radic  emitters.  He  wrote  the  program  in  extend¬ 
ed  FORTRAN  language  for  the  CEC  6600  digital  computer  sys¬ 
tem.  In  his  own  words: 

This  program  is  a  geometric  optics  model  of  wave 
propagation  through  an  inhomogeneous  atmosphere 
naving  a  vertically  stratified  index  of  rerrac- 
tion.  The  program  calculates  the  direction 
of  wavefront  propagation  by  solving  the 


Euler-La g r an ge  equations  cf  ra”s  normal  to  incre¬ 
mental  surfaces.  The  ray  trajectories  are  then 
used  to  coirpute  the  relative  emitter  field 
strength  cr  power  density  (normalized  to  tree 
space;  as  a  function  of  altitude  and  distance 
alone  the  earth's  surface.  Fields  which  are  re¬ 
flected  from  the  earth  are  attenuated  cy  a  Fresnel 
reflection  coefficient  and  a  surface  roughness 
factor.  The  elevation  ancle  and  time  of  propaga¬ 
tion  are  calculated  along  each  ray  patn  to  deter¬ 
mine  the  direction  of  the  wavefront  propagation 
vector  and  the  phase  relationship  between  inter¬ 
fering  wavefronts  for  the  field  strength  and  power 
density  computations.  [Ref.  39] 

Jim  Blake,  a  student  at  the  Naval  Postgraduate  School,  con¬ 
verted  this  program  for  operation  on  the  IBM  3933  computer 
system.  The  program  was  further  modified  for  use  in  this 
study  to  include  plotting  of  "M"  and  "d.M/dZ"  versus  height. 
The  model  is  applicable  to  propagation  above  30  MHz.  Given 
an  isotropic  emitter  cf  known  frequency,  polarization,  pulse 
width,  and  altitude  the  " N ",  "dN/dZ",  "M",  "dM/dZ"  and  free 
space  normalized  power  density  and  relative  field  strength 
are  calculated  as  a  function  cf  altitude  and  distance  along 
the  earth's  surface.  In  the  ray  trajectory  diagrams  an  ar¬ 
tificial  upward  curvature  of  the  rays  is  present  due  to 
plotting  the  earth's  surface  along  a  linear  rather  than 
curved  axis. 

The  following  parameters  were  utilized  in  running  the 
program  : 

Frequency  .  :  2900  MHz  or  9800  MHz 

Pulse  width .  :  t>.  5  ji sec  or  1340  ji sec 

Transmitter  height  (alove  sea  level):  90  ft 

(above  ground  level):  15  ft 

Emitter  Polarization  .  :  Vertical 

Earth  Surface  Type  .  :  Very  Dry  Land  (Type  2) 

Standard  Deviation  in  Height  .  :  Smooth  Plain  (Type  3) 

Figures  35,  39,  43,  47,  51,  55,  59  and  63  present  heignt 

gain  curves  computed  at  2900  MHz  for  a  15  foot  high  emitter 
(the  approximate  height  of  an  antenna  on  top  of  a  tracked 
vehicle).  Vertical  reference  axes  are  drawn  at  20  nautical 
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strength  relative  to  free  space  values.  A  scale  for  measur¬ 
ing  relative  fieli  strength  is  given  in  the  upper  right  cor¬ 
ner  of  each  plot. 

Wind  direction  £} ufology  (i.e.  N  for  north,  5  for  south, 
E  for  East,  etc.)  was  assigned  based  on  the  following: 

Radiosonde  Wind  Direction  Assigned  Direction 

Reading  (Degrees)  Symbol 


338 

-  22 

N 

22 

-  68 

N  E 

68 

-  1  12 

rp 

1  12 

-  153 

SE 

1  58 

-  2C2 

S 

202 

-  248 

SW 

00 

i  M 

-  252 

W 

282 

-  333 

NW 

B.  SURFACE  ANE  ELEVATED  DUCTS 

Atmospheric  •  Refractivity  provided  nine  (9)  distinct 
groupings  for  comparison.  The  first  eight  groupings  have 
some  form  of  ducting.  The  ninth  group  is  the  composite 
listing  of  all  radiosonde  launches  which  showed  no  ducting 
present.  The  three  sets  of  surface  based  ducts  caused  by  a 
surface  layer  are  defined  in  Figure  25.  The  numbers  depict¬ 
ed  are  the  averages  taken  from  the  data/findings  within  each 
group/set . 
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Set  1  is  pictoriaily  represented  tv  28  day  19  30302  in 
the  folio win  g  figures.  The  height  of  the  surface  based  duct 
averaged  approximately  1135  feet.  The  intensity  of  the  duct 
was  approximately  48  M-units.  There  were  too  few  data  re¬ 
cordings  to  determine  a  dominant  wind  direction  at  tae 
75-aCO  feet  and  300 0-4000  feet  levels.  A  Nor  t  i  i  w e  s  t  win  d  w a  s 
dominant  at  tne  43C0-5253  feet  level.  The  minimum  trapp^  i 
freyuency  for  this  set  was  51.3  NHr.  This  occurred  5r  of 
the  time  compared  to  the  15^  of  Figure  29. 

Figure  3u  shows  that  ray  trapping  occurs  mainly  between 
the  earth’s  surface  and  1130  feet.  The  area  above  the  duct 
shows  a  distinct  absence  of  rays.  This  would  be  consider  a  5 
the  radar  hole.  while  geometric  optics  predicts  that  there 
is  ao  field  present  i  r.  this  region,  geometric  optics  is  not 
able  to  solve  for  diffracted  fields  or  fields  resulting  from 
leak.y  modes  which  are  often  present  in  atmospheric  ducts. 
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Set  2  is  pictorial!/  represented  by  24  May  3  0  300  02  i  r. 
the  following  figures.  The  height  of  the  surface  based  duct 
averaged  approximately  393  feet.  The  width  of  the  duct  was 
approximately  22  M-units.  West  and  Northwest  winds  at  the 
75-400  feet  level.  North  and  Northwest  winds  at  the 
3000-4000  feet  level,  and  North  and  Northwest  winds  at  tne 
4000-5250  feet  level  were  dominant.  The  minimum  trapped 
freguency  for  this  set  was  80.2  MHz.  This  occurred  125  of 
the  time  compared  to  the  22%  of  Figure  28. 

Figure  38  shews  that  ray  trapping  occurs  mainly  between 
the  earth’s  surface  and  700  feet.  The  area  above  the  duct 
shows  a  radar  hole  extending  from  70  to  200  nautical  miles. 


is  pictorially  represented  by 


Set  3 


2  3  May  7  3  120  32  i 


the  following  figures.  The  height  of  the  surface  based  Juc  ■. 
averaged  approximately  214  feet.  The  width  of  tee  duct  was 
a  p  proxi  ma  tely  6  M-ur.its.  North  winds  at  the  75-400  feet 
level.  North  and  Northwest  winds  at  tne  1000-4010  feet  lev¬ 


el,  and  Northwest  winds  at  the  403  0-5250  feet  level  ert 
dominant.  The  minimum  trapped  frequency  for  this  set  was 
684.0  MHz.  This  occurred  18X  cf  tne  time  compared  to  the  2 r 


of  Figure  23. 

Figure  42  shew; 


that  ray  trapping  occurs 


the  earth's  surface  and  200  feet. 


?he  area  arove  the  duct 


shows  the  presence  of  a  radar  hole  extending  from  30  to  200 
nautical  miles. 


Set  4  is  yictorially  represented  by  26  May  3  0  000  02  in 

the  following  figures.  The  height  of  the  surface  Last!  hurt 
averaged  approximately  12^0  feet.  The  elevated  layer  was 
approximately  1093  feet  thick.  The  width  of  tne  duct  was 
approximately  49  M-units.  There  were  too  few  lata  record¬ 
ings  to  determine  a  dominant  wind  direction  at  tne  75-40  0 
feet,  3000-400C  feet  and  4  030-  5250  feet  levels.  The  minimum 


Set  5  is  pictorial!?  represented  t  y  _  '  i  7 
following  figures.  The  height  or  tu  surf  ,i_  -  •. 

eraged  approximately  cUO  feet.  T:.e  v  i:  : 

proximateiv  531  feet  thick.  The  writ'-.  >; 
pioximately  1  1  M-units.  There  were  t  ; o  :  s  •  1  . 

to  determine  a  dominant  wind  direction  a*  r  - 

3003-4000  feet  and  4030-5250  feet  lev-1. ; 
trapped  frequency  for  this  set  was  13  2.3 
2  5!  of  the  time  compared  to  the  15^  of  Figure 

Figure  50  shows  ray  trapping  ccc  jru;.  .  .  •.  • 

earth''s  surface  and  S00  feet.  There  is  ar.  ar  o 
arove  the  duct  (racar  hole)  extending  from  ’’J  *_  j  _ 
cal  miles  and,  a  low  altitude  hole  extending  fr  at 
nautical  miles.  The  hole  represents  the  to  gun  . 
earth's  horizon  where  rays  are  unable  to  pe  net  rat 
frequently  referred  to  as  the  earth's  shadow  region 
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Figure  U9:  ’Modified  Eefractivity  Profile 
2  May  7  8  00 00 Z  ' 


Set  6  is  pictoriaily  represented  by  17  Oct  7  9  0000Z  in 

the  following  figures.  The  height  of  the  surface  based  duct 
averaged  approximately  1330  feet.  The  elevated  layer  was 
approximately  879  feet  thick.  The  width  of  the  duct  was  a^ 
prcximately  26  K-units.  West  winds  at  tae  75-  40  3  feet  lev¬ 
el,  North  winds  at  the  300G-4QC0  feet  level,  and  North  and 
Northwest  winds  at  the  4000-5250  feet  level  were  dominant. 
The  minimum  trapped  frequency  for  this  set  was  45.7  MHz. 
This  occurred  13??  cf  the  time  compared  to  the  11%  of  Figure 
28. 

Figure  54  shows  ray  trapping  occurring  between  the 
earth’s  surface  and  1103  feet.  There  is  an  absence  of  rays 
above  the  duct  (radar  hole)  extending  from  60  to  230  nauti¬ 
cal  miles  and,  a  low  altitude  hole  extending  from  30  to  80 


nautical  miles 


Set  7  is  pictoriaily  represented  by  25  lay  SO  OOOOZ  in 
the  following  figures.  The  elevated  duct  averaged  approxi¬ 
mately  12S5  feet  in  thickness.  The  elevated  layer  was  ap¬ 
proximately  755  feet  thick.  The  width  of  the  iuct  was  ap¬ 
proximately  20  M-urits.  Tr.ere  were  too  few  data  recordings 
to  determine  a  dominant  wind  direction  at  the  75-400  feet, 
30CC-4000  feet  and  4000-5250  feet  levels. 

Figure  58  shews  the  presence  of  an  elevated  duct.  A  low 
altitude  (100  0  feet)  radar  hole  is  present  below  the  duct 
from  50  to  200  nautical  miles. 


Figure  56: 


’At ecs cheric  Refractivity 
25  May  3C  OOOOZ* 


Profile  - 
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’Modified  Refractivitv  Profile  - 
25  May  8  C  OOOOZ  ' 
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Figure  57 


Set  8  is  pictorially  represented  fcy  3  Oct  79  1200Z  i  r. 
tae  following  figures.  The  elevated  duct  averaged  approxi¬ 
mately  850  feet  in  thickness.  The  elevated  layer  was  ap¬ 
proximately  666  feet  thick.  The  width  of  tne  duct  was  ap¬ 
proximately  13  M-urits.  North  and  East  winds  at  tne  7  5-4  J  O 
feet  level,  North  winds  at  the  3000-4000  feet  xevel,  ar.i 
North  and  Northwest  winds  at  the  4000-5250  feet  level  were 
do  minan t . 

Figure  62  shews  the  presence  of  an  elevated  duct.  A  low 
altitude  (2000  feet)  radar  hcle  is  present  below  the  duct 
from  50  tc  200  nautical  miles. 


Figures  64  thru  69  attempt  to  show  the  correlation  of 
wind  direction  and  wind  speed  tc  the  estabiisnment,  location 
and  heights  of  ducts.  Unfortunately,  wind  direction  and 
wind  speed  data  were  very  limited  between  400  feet  and  3000 
feet.  No  distinctive  pattern  was  able  to  be  obtained  from 
the  results.  To  compare  the  wind  speeds  depicted  here  with 
Figure  20  multiply  meters  per  second  by  two  to  get  knots 
(1  m/s  =  1.945  knots). 
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Surface  based  ducts  caused  by  a  surface  layer  occurred 
35%  of  the  time;  surface  based  ducts  caused  r-y  an  elevated 
layer  occurred  16%  of  the  time;  elevated  ducts  occurred  13% 
of  the  time;  and,  no  ducts  were  present  39%  of  the  time. 
During  May  at  OOOOZ,  surface  based  ducts  caused  oy  a  surface 
layer  were  predcrinant  occurring  52%  of  the  time.  Of  this 
52%,  Group  2  occurred  58%  of  the  time.  During  May  at  12092, 
there  were  no  signs  cf  ducting  52%  of  the  time.  However 
surface  based  ducts  caused  by  a  surface  layer  did  occur  37% 
of  the  time.  Cf  this  37%,  Group  3  was  present  68%  of  tne 
time.  During  Cctober  at  OOOOZ,  surface  based  ducts  caused 
by  elevated  layers  occurred  38%  of  the  time.  Of  this  381, 
Group  6  occurred  £5%  of  the  time.  During  October  at  12032, 
there  was  no  sign  of  ducting  56%  of  the  time.  Again  how¬ 
ever,  surface  based  ducts  caused  by  a  surface  layer  aid  oc¬ 
cur  24%  of  the  time  of  which  86%  were  Group  3. 


C.  LAND-SEA  BE  EE ZE 

Thru  analysis  it  appears  that  between  the  heights  of 
951-1200  meters  is  where  the  weather  maps  used  in  Chapter 
III  get  there  readings  for  wind  direction  and  speed  in  the 
Dhahran  area.  Figures  70,  71,  and  72  contain  the  data  re¬ 
corded  during  the  radiosonde  launches.  This  information 
matches  pretty  well  with  Figures  18,  19  and  20. 
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Figure  70;  'May  1978,  1  979,  t  1  930  Wind  Direction  (%) 

at  3120  -  3937  feet' 
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Figure  71:  'October  1978,1  979,  5  1930  Wind  Direction 

(K)  at  3120  -  3937  feet' 
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However,  this  data  dees  not  shew  the  effect  of  the  land-sea 
breeze.  Data  between  the  heights  of  23  to  122  meters 
(75-400  feet)  shows  a  better  picture  (See  Figures  41,42  & 
43). 
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Figure  73: 


'Hav  1  978,  1  979,  £  1  930  Wind  Direction  {%) 
*  at  75  -  4 CC  feet • 
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Figure  75:  'May  and  October  1978,1979.  o  1930  Kind 
Speed  (m/sec)  Averages  at  75  -  400  feet' 
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At  000 0Z  winds  frcm  the  west  at  2.7  meters  per  second 
(5.4  knots)  form  the  land  breeze.  At  1200Z  winds  from  the 
north  at  7  to  S  meters  per  section  or  from  the  east  at  4  me¬ 
ters  per  second  (8  knots)  form  the  sea  breeze.  When  the 
wands  blow  from  the  west,  northwest  or  southwest  forming  tne 
land  breeze,  they  shifted  to  the  north,  northeast  or  east  by 
1203Z  to  form  the  sea  breeze.  Wind  speed  was  stronger  for 
the  sea  breeze  as  expected.  No  correlation  could  be  deter¬ 
mined  between  the  occurrence  of  land-sea  breezes  and 
d  ucting . 


V.  2 AC TI CAL  APPLICA  TIC  NS 

Tactical  military  commanders  car  not  survive  without  mak¬ 
ing  use  of  and  exploiting  the  EM  spectrum.  Although  the  EM 
spectrum  has  been  utilized  the  military  has  not  fully  con¬ 
sidered  atmospheric  effects  on  the  EM  spectrum.  The  follow¬ 
ing  commonly  used  systems  exemplify  the  affects  atmospheric 
anomalies  have  or.  EM  propagation. 

A.  RADAR 

The  military  has  been  the  major  user  of  radar  and  the 
contributor  of  its  developmental  cost.  The  major  areas  of 
radar  application  for  the  military  includes  but  is  not  lim¬ 
ited  to  air  traffic  control,  aircraft  navigation,  snip  safe¬ 
ty,  remote  sensing  (i.e.  used  as  a  remote  sensor  of  the 
weather  or  as  an  ionospheric  sounder)  ,  surveillance  and  for 
control  and  guidance  cf  weapons.  Conventional  radars  genec- 
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The  radar  equation  relates  the  range  of  a  radar  to  the  char¬ 
acteristics  of  the  transmitter,  receiver,  antenna,  target 
and  environment . 

~Pr  G_  Ae  cr  I-  (n)  “  1/4 

=  - -  (^n  10) 

K  (T„  ♦  Te  )  B  0  L  ^  (S/N),^_ 

where,  RWAv  =  maximum  radar  range  (m) . 

Pr  =  peak  power  output  of  radar  (watts). 

Gr  =  antenna  gain  cf  radar. 

Ad  =  effective  antenna  area  (m^)  . 

cr  =  radar  cross  section  (m2)  . 

I L  (n)  =  integration  improvement  factor. 

K  =  Boltzmann's  constant  (1.38  x  10'23  J/Deg) 

Tt_  =  antenna  noise  temperature  (°K). 

Te  =  equipment  noise  temperature  (“K). 

3n  =  noise  bandwidth  (Hz)  . 

=  system  losses. 

(S/N)i(^  =  minimum  signal  to  noise  ratio  of  a  single 
pulse. 

I  his  equation  finds  the  maximum  radar  range.  However  it 
does  not  take  into  consideration  atmospheric  refractivity 
which  as  seen  in  Figures  34,  38,  42,  46,  50,  54,  58  and  52 
can  have  a  severe  impact  on  what  a  radar  actually  sees. 

Battlefield  surveillance  radars  such  as  tne  AN/PPS-5 
which  can  detect  men  cut  to  50C0  meters  and  vehicles  out  to 
10,000  meters  or  the  AN/PPS-6  which  can  detect  personnel  out 
to  1500  meters  and  vehicles  out  to  3000  meters  ioox  out  hor¬ 
izontally  along  the  earth's  surface  and  therefore  will  not 
he  affected  by  atmospheric  refractivity.  Aircraft  control, 
precision  approach  (AN/TPN-25  S  AN/TP N-13A)  and  short  range 
air  defense  (AN/VPS-2)  radars  which  look  out  to  a  range  of 
40  miles  will  te  slightly  affected.  Atmospneric  conditions 
were  present  6.45!  cf  the  time  which  would  cause  EM  waves  to 


bend  (see  Figures  34  r  46)  creating  a  radar  hole  above  the 
2000  foot  level.  Long  range  surveillance  (A'J/TPS-32)  and 
air  defense  (AH/TPS-43E)  radars  which  have  ranges  in  excess 
of  240  nautical  miles  (r.mi)  will  be  severely  affected.  Ra¬ 
dar  hoi ^s  were  present  starting  at  a  distance  of  30  nai  from 
the  radar  at  a  height  of  200C  feet  57%  of  tn e  time  under 
s  t  udy . 


E.  COM  HUSICATICNS 

H?  and  VHF  communication  systems  operating  at  frequencies 
less  than  45.7  MHz  would  not  have  been  appreciably  affected 
by  atmospheric  refractivity  during  the  study  months.  Trans¬ 
mitters  operating  above  the  following  frequencies  would  have 
experienced  extended  ranges  at  the  also  listed  percentages: 


Frequency 
45.7  MHz 

47.3  MHz 

51.3  MHz 
80.2  MHz 

132.3  MHz 
684. C  MHz 


Extended  Range  Percentage 
13% 

15% 

20% 

32% 

34% 

52% 


The  trapping  cf  EM  signals  in  surface  ducts  and  elevated 
ducts  alsc  showed  that  signal  strength  loss  was  not  as  se¬ 
vere  as  it  would  have  been  had  there  been  no  duct  present. 
Hence  the  intercept  of  signals  could  be  made  at  longer  dis¬ 
tances  and  with  less  sensitive  receivers. 
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VI.  CC NCI OS I OSS  AND  RECOMMENDATIONS 

Surface  based  ducts  and  elevated  ducts  were  present  6 1  f. 
of  the  time  during  the  months  of  May  and  October  in  1973, 
1979  and  1960.  These  ducts  wculd  most  severely  affect  long 
range  surveillance  radars  due  to  the  presence  of  radar 
holes,  and  communication  systems  by  extending  their  normal 
transmitting  range.  Radar  holes  were  predominantly  present 
starting  at  a  distance  of  SO  nmi  from  the  radar  at  a  height 
of  2000  feet  ar.d  rising.  Elevated  d  ucts  present  10%  of  the 
time  to  a  maximum  height  of  2000  feet  will  also  degrade  the 
performance  of  Side  looking  Airborne  Radar  (SLAR)  and  Joint 
Army  Air  Force  Surveillance  and  Attack  Radar  System  (JSTAF.3) 
(aircraft  height  approximately  15,000  feet)  if  they  were  to 
be  utilized  in  this  area  of  the  world. 

The  land-sea  breeze  phenomenon  was  present  during  the 
study  months.  Nc  correlation  could  be  determined  between 
the  occurrence  cf  the  land-sea  breezes  and  ducting. 

A  portion  of  this  thesis  attempted  to  show  the  correla¬ 
tion  of  wind  direction  and  wind  speed  to  tne  establishment, 
location  ar.d  heights  of  ducts.  Wind  direction  and  speed 
data  was  very  limited  between  400  feet  and  3000  feet.  By 
coincidence  this  area  happened  to  be  located  eitner  directly 
inside  or  just  above  the  duct.  Hence  no  distinctive  pattern 
was  able  to  be  obtained  from  the  results.  In  the  future, 
more  data  is  needed  to  be  obtained  from  radiosonde  launches 
in  the  predicted  duct  area  in  order  for  a  more  reasonable 
analysis  to  occur. 
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APPENDIX  A 


COMPOTES  PROGRAM 

//PETERSEN  JOB  ( 2  3  10 , 0  39  1 )  ,  '  SM  C  23  10 •  , CLAS5=C 
//♦MAIN  CF.G=NPGVM  1 . 2 3 1 0? , L  IN ES  =  (60) 

//♦FORMAT  PR,  D  P  NAM  E  =  PLOT .  SYS7ECTR/DEST=I0CAL 
//  EXEC  FF.TXC  I G  E  ,  PA  F  M.  FOF.T=  ’  MAP,  GOST  MT, XBE r * 

// FORT. SYS IN  DC  ♦ 

C  THIS  PROGRAM  WAS  ORI3INAIIY  DEVELOPED  BY  RAYMOND  P.  .. 

C  WASKY.  PROGRAM  CONVERTED  FOR  USE  ON  THE  IBM  3033  BY.  . 

C  JIM  EL  A  E  E .  PROGRAM  MODIFIED  FOR  "M"  AND  "DM/DZ"  EY  .. 

C  3  Y  WAYNE  F .  PETERSEN. 

C - THIS  PROGRAM  IS  A  GEOMETRIC  OPTICS  MODEL  OF  WAVE 

C - PROPAGATION  THROUGH  AN  I NKOMCGE NSO OS  ATMOSPHERE  HAVING 

c - A  V  ERTIC  AIL  Y  STRATIFIED  INDEX  OF  REFRACTION.  THE 

C - PROGRAM  CALCULATES  THE  DIRECTION  CF  WAVEFRONT 

C - PROPAGATION  EY  SOLVING  THE  EULE R-LAGRANGE  EQUATIONS  OF 

C - RAYS  NORMAL  TO  INCREMENTAL  WAVEFRONT  SURFACES.  THE  .. 

C - RAY  TRAJECTORIES  ARE  THEN  USED  TO  COMPUTE  THE  RELATIVE 

C - EMITTER  FIELD  STRENGTH  OR  POWER  DENSITY  (NORMALIZED  IC 

C - TO  FEES  SPACE)  AS  A  FUNCTION  DF  ALTITUDE  AND  DISTANCE 

C - ALONG  THE  EAFIH’S  SURFACE.  FIELDS  WHICH  ARE  REFLECTED 

C - FROM  THE  EARTH  ARE  ATTENUATED  BY  A  FRESNEL  REFLECTION 

C - COEFFICIENT  AND  A  SURFACE  ROUGHNESS  FACTOR.  THE 

C - ELEVATION  ANGLE  AND  TIME  CF  PROPAGATION  ARE  CALCULATED 

C - ALONG  EACH  RAY  PATH  TO  DETERMINE  THE  DIRECTION  OF  THE 

C - WAVEFRONT  PROPAGATION  VECTOR  AND  THE  PHASE  E  EL  A  TIC  N- 

C - SHIP  BETWEEN  INTERFERING  WAVEFRONTS  Fun  THE  FIELD 

C - STRENGTH  AND  POWER  DENSITY  COMPUTATIONS. 

C 
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C  DATA  INPUT  * 

C  * 

C  THIS  PROGRAM  WILL  ACCEPT  TWO  TYPES  OF  DATA  INPUT.  SET  * 

C  DAICTF.L  =  3  FCE  INPUT  OF  HEIGHT  AND  RSFRACTIVITY .  SET  * 

C  SET  DAT CTL  =  1  FCE  INPUT  OF  HEIGHT,  N,  M ,  DNDS,  l  DMDZ.  * 

C  INPUT  FCE  MAT  IS  6(F10.  3).  IF  THE  HIGHEST  HEIGHT  IN  FEET  * 

C  IS  LESS  THAN  YCUE  ASSIGNED  Y  MAX  you  can  only  request  * 

C  PROFILES  AND  PLOTS  FOE  "N",  "IN/DZ",  THE  ray  trace  ar.  a  * 

C  POWER  DENSITY.  SEE  SUBROUTINE  REFECT  FOR.  explanation  * 

C  * 

C  CHECK  READ  STATEMENTS  FOR.  OTHER  NECESSARY  DATA  TO  EE  * 

c  input.  * 

C  * 


c 


C  PLOTS 

C  _ 

C  THERE  ARE  SIX  PLOTS  AVAILABLE  IN  THIS  PROGRAM  : 
C  1.  R  EFRACT IV ITY  (N- UN  ITS)  VS  HEIGHT  (FEET) 

C  2.  DN/DZ  { N-U  NIT  S/KM)  VS  HEIGHT  (FEET) 

C  3.  MODIFIED  INDEX  OF  REFFACTIVITY  (M-UNIIS) 

c  vs  Height  (feet) 

C  4.  DM/DZ  (M- UNIT S/KM)  VS  HEIGHT  (FEET) 

C  5.  RAY  TRACE 

C  6.  RELATIVE  POWER  DENSITY 

C 


* 

* 

* 

* 

$ 

* 

* 

* 

* 

* 

* 


£*$$***  *##$*#*$***$***:********:*******  ***$*$*#***#***$:$**#**  * 

C 

c 
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*  * 


c 

c 

C  NCASE 
C  :I  PRO 
C 
C 

C  NDAT  A 
C 

C  K  N  ( 1 ) 

C 

C  HN(NDATA) 
C 

C  HF  (? ) 

C 

C  RN  ( 1  ) 

C  RN  (N  DATA) 
C  AHS 
C 

C  AHO 
C 

C  YMIN 
C 

C  YMAX 
C 

C  XDELTA 
C 

C  XFINAL 
C 

C  ELOS1 
C  ELOS2 

C 


VARIABLE  DEFINITIONS  * 

* 

-  NUMBER  OF  FROG RAM  RUNS.  * 

-  CODE  FOR  R EFRACTI VII Y  MODEL  * 

1=0  FREE  SPACE  MODEL;  * 

=  1  EXPONENTIAL  MODEL;  =2  INPUT  MODEL)  * 


-  NUMBER  OF  REFRACT  IV  IT Y  PROFILE  DATA  * 

LEVELS  (NEATA  =  1  IF  NPRC  =  0  OR  1 )  * 

-  HIGHEST  ALTITUDE  IN  REFRACT  IV  IT Y  * 

PROFILE  (FEET)  * 

-  LOWEST  ALTITUDE  IN  F.  EF  RAC  II  VI T  Y  * 

FROFILE  (FEET)  * 

-  ALTITUDE  IN  REFRACTIVITY  PROFILE  * 

(MEIERS)  * 

-  REFRACTIVITY  AT  HN(1)  * 

-  REFRACTIVITY  AT  HN(NDATA)  * 

-  EARTH  SURFACE  ALTITUDE  IN  FEET  ABOVE  * 

-  SEA  LEVEL  * 

-  EMITTER  ALTITUDE  IN  FEET  ABOVE  * 

-  SEA  LEVEL  * 

-  MINIMUM  ALTITUDE  IN  FEET  ABOVE  SEA  * 

LEVEL  FOR  PRINT  AND  PLOT  OUTPUT  * 

-  MAXIMUM  ALTITUDE  IN  FEET  ABOVE  SEA  * 

LEVEL  FOR  PRINT  AND  PLOT  OUTPUT  * 

-  DISTANCE  INTERVAL  IN  NAUTICAL  MILES  * 

FOR  PRINT  AND  PLOT  OUTPUT  * 

-  MAXIMUM  DISTANCE  IN  NAUTICAL  MILES  * 

FOR  PRINT  AND  PLOT  OUTPUT  * 

-  HIGHEST  RAY  ANGLE  IN  DEGREES  * 

-  LOWEST  RAY  ANGLE  IN  DEGREES  * 

(DIFFERENCE  BETWEEN  ELOS  1  AND  ELOS2  * 
CAN  BE  NO  MORE  THAN  ONE  DEGREE)  * 

-  EMITTER  FREQUENCY  IN  MEGAHERTZ  * 


8  1 


c 

PW 

-  EMITTER  PULSE 

WIDTH  IN  MICROS  EC 

* 

c 

(IF  EMITTEF  IS  CONTINUOUS  WAVE  SET 

$ 

c 

PF  =  10  00  COO. 

0  ) 

* 

c 

N  HV 

-  CCDE  FOR  EMITTER  POLARIZATION 

* 

c 

(=1  HCRIZ;  =2 

VERT) 

* 

c 

NSL 

-  CCDE  FOR  EARTH 

SURFACE  TYPE 

* 

c 

=1  SEA  WATER; 

=2  VERY  DRY  LAND; 

* 

c 

=  3  AVERAGE  LAND;  =4  VERY  MOIST  LAND) 

* 

c 

N  SMS 

-  CCDE  FOR  SURFACE  ROUGHNESS 

* 

c 

(SEE  TABLE  BE 

LOW) 

* 

c 

CODE 

STANDARD  DEVIATIONS  OF  HEIGHT 

* 

c 

(NF.MS) 

(SEA) 

(LA  ND) 

* 

c 

0 

o 

• 

o 

0 

* 

c 

1 

CN 

• 

o 

9 

* 

c 

2 

0.6 

30 

* 

c 

3 

1.  1 

56 

* 

c 

4 

1.7 

1  12 

* 

c 

5 

2.6 

2  14 

* 

c 

6 

4.  3 

42S 

* 

c 

7 

8.6 

12  88 

* 

c 

8 

12.  9 

2146 

* 

KB2F 

KGRAD 

KRAY 

KPLOT 


KM  IB 


-  CCDE  FOR  REFRACTIVITY  PROFILE  PRINTOUT  * 

(=0  NO  PRINTOUT;  =1  PRINTOUT)  * 

-  CCDE  FOR  RFFRACTIVITY  GRADIENT  PROFILE  * 

PRINTOUT  (=0  NO  PRINTOUT;  =1  PRINTOUT)* 

-  VARIABLE  NOT  USED  -  SET  =  0  * 

-  CCDE  FOR  RELATIVE  FIELD  STRENGTH  OR  * 

POWER  DENSITY  PRINTOUT  * 

(=0  NO  PRINTOUT;  * 

=  1  RELATIVE  FIELD  STRENGTH  PRINTOUT;  * 
=2  RELATIVE  POWER  DENSITY  PRINTOUT)  * 

-  CCDE  FOE  MODIFIED  INDEX  OF  REFR A CTI VII Y  * 

PRINTOUT  !=0  NO  PRINTOUT;  =1  PRINTOUT)* 
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KMGRAD  -  CODE  FOR  MCE  INDEX  OF  RZFRAC  GRADIENT  * 

FRINTOUT  (=0  NO  PRINTOUT;  =1  PRINTOUT)  * 
N  REF  -  CODE  FOR  R  EFRACTI VIT  Y  PLOT  * 

(=0  NO  PLOT;  =1  PLOT)  * 

NGEAD  -  CODE  FOR  R  EFRACTI VIT Y  GRADIENT  PLOT  * 

(=0  NO  PLOT;  =1  PLOT)  * 

N  F  A  Y  -  CCDE  FOR  RAY  TRACE  FLOT  * 

(=0  NO  PLOT;  =1  PLOT)  * 

N  PL  j  T  -  CCDE  FOE  RELATIVE  FIELD  STRENGTH  OR  * 

POWER  DENSITY  PLOT  * 

(=0  Nc  plot;  * 

=  1  RELATIVE  FIELD  STRENGTH  PRINTOUT;  * 
=2  RELATIVE  POWER  DENSITY  PLOT)  * 

N M I R  -  CCDE  FOR  MODIFIED  INDEX  OF  HZFRACTIVITY * 

D LOT  (=0  NC  PLOT;  =1  PLOT)  * 

NMGR AD  -  CCDE  FOR  MCE  INDEX  OF  REFRAC  GRADIENT  * 

FLOT  (=0  NC  PLOT;  =1  PLOT)  * 

SCALE  -  SCALE  FACTCP  FOR  ENLARGING  OR  REDUCING  * 

FLOT  SIZE  FROM  ITS  NORMAL  5  X  10  * 

INCH  FORMAT.  NO  R  MAI  PLOT  SIZE  GIVEN  * 

WITH  SCALE  =1.0  * 

IAHD  -  NUMBER  OF  CHARACTERS  IN  PLOT  BANNER  * 

A  HD  -  CHARACTERS  IN  PLO  I  BANNER  * 

I XTL  -  NUMBER  OF  CHARACTERS  IN  X-AXIS  LABEL  * 

XTL  -  CHARACTERS  IN  X-AXIS  LABEL  * 

I YTL  -  NUMBER  OF  CHARACTERS  IN  Y-AXIS  LABEL  * 

YTL  -  CHARACTERS  IN  Y-AXIS  LABEL  * 

ITTL  -  NUMBER  OF  CHARACTERS  IN  FIRST  LINE  OF  * 

PLOT  TITLE  * 

TTL  -  CHARACTERS  IN  FIRST  LINE  OF  PLOT  TITLE  * 

ITLE  -  NUMEER  OF  CHARACTERS  IN  SECOND  LINE  OF  * 

tLOT  TITLI  * 

TLE  -  CHARACTERS  IN  SECOND  LINE  OF  PLOT  TITLE* 


c 

IMPLICIT  BEAI*8  (A-H,0-Z) 

CCMMCN  /CNECCM/  CRH, CRX , CEG, DTR , PH AS E,T HSTA , C, C  1 , C2 

1  ,RHO  (51,30)  ,PHI  (51  ,  3C)  ,?I,  CMF,CNAUT, ienot 

2  ,FR£,TPS, AESRH, NHV, KSL, NRMS,JAEEA, NDATA 

CCMMCN  /TWOCCM/  EA  ,  DELX  ,  AHS ,  AHO ,  ELO  ( 5  1)  , E  XMAX,  XFI NAI  . 

1  ,jray 

COMMON  /THE  CCM/  VX  (3  0)  ,  Y  MIN,  YM  A  X,  XDI V  ,  YDI  V 

1  ,  H  (51,30)  ,0(51,30)  ,  HN  (100)  /EN  (10  0)  ,  S  (51/30) 

2  , JPLT,KEEF  ,K3RAE, KP1CT, JPLCT, IPLACE, JCASE, NZLO 

3  ,NREF, NBA  Y,NPRO,NPLCT,NGF.  AD, SCALE 

CCMMCN  /FCECCM/  DNDZ  (1  00  )  ,RM  (10  0)  ,  DMDZ  ( 10  0)  ,  KMIE 
1  ,  kmgrad,NMIR,NMGRAE,HF  (100)  ,DATCTL,  Z,LL 

INTEGER  DATCII 

EEAI  +  16  PCL  (2)  ,  TER  (4)  ,  REEMOD  (2,  3) 

DATA  POL/  10  HRCRIZONT  AL, 1 0 EVERT I  CAL  / 

DATA  TEE/  10  ESEA  WATER  , 1 OEDRY  GROU ND , 10  HA VG  GROUND, 

1  10KKET  GROUND/ 

DATA  REFMCD/ 10HFP.EE  SPACE,  10H  MODEL  ,  1  0 HEXPO NENTIA  , 

1  10HL  MODEL  ,10HINPUT  PROF,10HILE  MODEL  /• 

DATA  CM/2.998D+08/ 

C  SET  UP  THE  NUMBER  OF  PROGFAM  RUNS 

R EAE  300 , NCASE 

C  IF  LATA  CCNTECL  (DATCTL)  =  1,  THEN  INPUT  DATA 

C  CONTAINS  HEIGHT,  N,  M,  DNEZ,  S  DMDZ.  IF  =  0,  THEN  .. 

C  INPUT  DATA  CONTAINS  HEIGET  5  N  ONLY. 

READ  300, EATCTL 
DO  200  IC AS E  =  1 , NCASE 
C 

C  READ  AND  ERINT  THE  INPUT  LATA 

R  EAE  300, NPRC 
R  EAE  300  , NDATA 


c 

c 


DO  £  ID=1, NEATA 
IF  (CATC  T I)  4,4,3 

3  READ  320,  Hr  (ID)  ,EN  (ID)  ,  R  H  (ID)  ,D  NDZ  (ID)  ,DMDZ  (ID) 

30  TC  5 

4  SEAT  320,  FF  (ID)  ,EN  (ID) 

CONVERT  HEIGHT  CATA  FROM  METERS  TO  FEET 

5  H N  (IE)  =  HF(ID)  *  3-  280840 
8  C  CNTIN'J  E 

REAL  32  0  ,  A  E5 ,AHC 
REAE  320,YMIN,Y MAX 
READ  320  ,XDEITA,XFINAL 
READ  320,EIOS1,  ELOS2 
REAE  320,FR;,PW 
REAE  300, NHV,NSL,NRMS 

REAE  300,KEEF,KGRAD,KRAY,  KPLOT,  KMI3,  KMGRA  D 
REAE  30 0, NREF, NGRA  D, NR  AY  ,  NPLOT,  NMIR, NMGRAD 

PRINT  THE  INPUT  DATA 
PRINT  330 

PRINT  340,NCASE,NPRO,  ( RE  FMOD  (IR  M,  N  P  RO  + 1 )  ,  IRM=1,2)  , 
1  NEATA 

PRINT  3  5  0  ,  A  HS,  A  HO,  YM IN  ,  Y  MAX 
PRINT  360 ,XCEITA,XFINAL 
PRINT  3  7  0  ,  ELCS  1 ,SLOS2 

PRINT  3  30,  FR£,  PW,POL  (NHV)  ,TSR  (NSL)  ,  NR  MS 

PRINT  390 ,KREF, KGRAD, KRAY, KPLOT 

PRINT  3 9 5 , KHIR , KMG R AD 

PRINT  4  0C,NREF,NGF.  AD,N  RAY,  NPLOT 

PRINT  405, NMIR,  KilGRAE 

STILL  CHECKING  INPUTTED  EATA 

PRINT  500 

DO  13  I  D  =  1,  NEATA 

IF  (EATCTI)  12,12,11 
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11  PRIST  51C,HF  (ID),HN{ID)  ,RN(ID)  ,  EM  (ID)  ,DNDZ(ID)  , 


C 

C 


C 

C 


C 

C 


c 

c 


i  cmdz(id) 

GO  IC  13 

12  PRINT  530, HF  (ID)  ,HN  (ID)  ,  BN  (ID) 

13  CONTINUE 

SET  UP  INITIAL  ATMOSPHERIC  REFR ACT IV IT  Y  CONSTANT 
IF  (NPRO-1)  14,15,20 

FREE  SPACE  MODEL 

14  C1=0.0 
C2=  C.  0 
GO  TC  25 

EXPONENTIAL  MODEL 

15  C 1 =3 13. 0 
C2=0. 00004386 
GO  TC  25 

PIECE-WISE  LINEAR  MODEL 
20  C2=0. 00004386 

C  1  =  RN  (1)  *DEXF  (C 2*HN  (  1}  ) 

SET  UP  INITIAL  CONDITIONS 
25  C=CMF*CM 
N EL C=50 

DENCT=1  ./EIOAT  (NELO) 

DELX=CNAUT*XDELTA/10. 

E  XM  AX=C  NAUT*XFINAL 
FRQ=1000000.C*FEQ 
P  LAC E=EX ft AX/DE IX 
TPW=PW/1 000CC0. 0 
ICAIC=K?ICT*NRA Y+NPLCT 
IPLACE=IFIX  (SNGL  (PLACE)  )  +1 
I PLCT=NREF+NGRAE+NR AY+NPICT 
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JCA  £E=IC  A  £E 
J  ?L 1=0 
N  ELC=NELC  + 1 


CALI  EOCJT  INES  I F  THERE  IS  A  PRINTOUT  OR  PLOT  Or  TH 
REFEACTIVITY  PROFILE 
J  PL  CT= 1 

IF  ( NREF )  30,30,31 

30  CONTINUE 

IF  (KR2F)  33,33,32 

31  CALI  PLOTTF 

32  CALI  REFECI 

CALL  ROUTINES  IF  THERE  IS  A  PRINTOUT  OR  PLOT  OF  TH 
REFEACTIVITY  GRADIENT 

33  J  PL CT=2 

IF  (NGRAD)  34,34,35 

34  CONTINUE 

IF  (KGRAE)  37,37,36 

35  CALI  PLCTTE 

36  CALI  REFECT 

CALL  ROUTINES  IF  THERE  IS  A  PRINTOUT  CP.  PLOT  OF  TH 
MODIFIED  INDEX  OF  KEFRACIIVITY  PROFILE. 

37  JFLCT  =  5 

IF  (NMI R )  50,50,51 

50  CONTINUE 

IF  (KMIR)  53,53,52 

51  CALI  PLOTTE 

52  CALL  REFECT 

CALI  ROUTINES  IF  THERE  IS  A  PRINTOUT  OR  PLOT  OF  TH 
MODIFIED  INDEX  OF  REFEACTIVITY  GRADIENT. 


53  J  PLOT 
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IF  (NMGRAD)  54,  54,  55 

54  CONTINUE 

IF  (KMGBAE)  57,57,56 

55  CALI  PLOTTF 

56  CALI  REFECT 

r* 

C  CHECK  IF  THZFE  ARE  ANY  FURTHER  CALCULATIONS 

57  IF  (ICALC)  200,200,40 
C 

C  CHECK  IF  THERE  IS  TO  BE  A  PLOT  OF  THE  RAY  TRACES. 

C  IF  SC,  CAIL  ROUTINE  TO  SET  rJP  THE  PLOT  AXES. 

40  IF  (NRAY)  44,44,42 
42  J  PL  CT=3 

CALI  PLOTTE 
44  CONTINUE 
C 

C  SET  U?  A  LOOP  TO  CALCULATE  THE  ALTITUDE  PROFILE  OF 

C  EACH  RAY. 

DO  100  1=  1 , NELC 
C 

C  INITIALIZE  FCR  INTEGRATION 

CRH=AKO 
CRX  =  0.  0 
ICCC=I-  1 

E  LO  (I)  =ELCS1  -DE NOT* FLOAT  (ICCC) 

EANG  =  DTR*EIC  (I) 

CRG=  ( (RA  +  CFH)/EA)  *DSIN  (BANG)  /DC  OS  (  EANG) 

JR A  Y  =  I 
C 

C  CALI  ROUTINE  TO  COMPUTE  RAY  TRACES  AND  PROPAGATION 

C  TIMES. 

CALI  RKINTG 
100  CONTINUE 
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RESET  PARAMETER  VALUES 
D  EL  X  =  1 0 . *DELX 
I?LACE=IPLACE/10+1 

CALL  ROUTINES  IE  THERE  IS  A  PRINTOUT  OR  PLOT  OF  TH 
RELATIVE  II  ELL  STRENGTH  OR  RELATIVE  POPER  DENSITY  . 
120  J  ?L CT=4 

IF  (NPLOT)  13  C  ,  130,  140 
130  CONTINUE 

IF  (KPLOT)  200 , 200,  1  50 
140  CALL  PLOTTE 
150  CALL  HGAIN 
200  CONTINUE 

CALL  LI3RARY  ROUTINE  FOR  CN-LINE  PLOTTING 
IF  (IPLOT)  220,220,2  10 
210  CALI  PLOT  (0.,0.,999) 

PRINT  430 
220  CONTINUE 


300  FCRKAT 

320  FORMAT 

330  FORMAT 

(615) 

(  6  F  1  0 

( 1  HI  , 

.3) 

2X,  'ATMOSPHERIC  RADIO  REFRACTIVE 

1 

'  COMPUTATION  S' 

//) 

340  FORMAT 

(2X,  ' 

NCA  S  E  =  * 

#16/ 

1 

2  X , ' 

NPRO  =  ' 

#16/ 

2 

2X,  » 

MOD  EL  =  » 

, 5X,2A10/ 

3 

2 X,  ' 

NDAT  A  =  ' 

, 16, 4X, 'LEVELS') 

350  FORMAT 

( 2  X,  ' 

A  H  S  =' 

,F10.2,4X,'FEET'/ 

1 

2 X,  • 

A  HO  =  ' 

,F1C.2,4X,'FEET'/ 

2 

2  X , ' 

YMI  N  =  ' 

,F10.2,4X,'FEET'/ 

3 

2 X,  ' 

YM  A  X  =» 

,F1C.2,4X, 'FEET') 

360  FORMAT 

( 2  X ,  ' 

XEELT A=' 

,  F  10.  2,4X  ,  '  NAUT  MI'/ 

1 

2 X,  ' 

XFINAL=' 

, F 10. 2,4X , ' NAUT  MI') 

370  FORMAT 

(2X,’ 

ELO S  1  ='  ,F10.2,4X,  'DEG'/ 

1 

2X,  • 

ELO  S2  ='  , F  10. 2,4X , ' DEG ' ) 

330  FORMAT 

(2X,  • 

FRQ  =  '  ,F  10.  2,4X  ,  *MHZ  •/ 

1 

2X,» 

PK=’ , F  10. 2,4X ,  ' MICROS  EC ’/ 

2 

2X,  « 

POL  AR  =  '  ,  5  X,  A 1  0/ 

3 

2  X ,  ' 

TERRAINS  , 5X,A1Q/ 

4 

2 X/  ' 

NRM5  =  *  ,16) 

390  FORMAT 

(  2  X  /  • 

KR  E  F  =  l  ,16/ 

1 

2  X ,  * 

KGR  AD  =  f  ,  16/ 

m 

2 

2  X ,  ' 

KR  A  Y  =  '  ,16/ 

m 

3 

2X,  » 

KPL  OT  =  '  ,16) 

• 

395  FORMAT 

{ 2  X  /  » 

KMIR  =  »  ,16/ 

1 

2X,  ’ 

KKGR  AD='  ,  16) 

400  FORMAT 

( 2X,  • 

NR  E  F  = '  ,16/ 

1 

2 X,  • 

NGR AD  =  *  ,  16/ 

2 

2X,  » 

NR  AY  =  *  ,16/ 

3 

2  X ,  ’ 

N?LOT  =  *  ,16) 

405  FORMAT 

( 2  X  ,  * 

NMIR=*  ,16/ 

i  *  < 

1 

2  X,  1 

N  MGR  AD  =  '  ,  16) 

430  FORMAT 

( 5  X,  1 

END  OF  FILE  CN  PLOTTER  TAPE') 

• 

500  FOR  M  AT (/ 3 X, ' 

HT  (M)  '  ,  4X  ,  •  HT  (FT)  *  ,  9  X,  *  N  1  ,  1  0  X,  '  M '  ,  6  X, 

• 

1  •  DNE2’,6X,*DMDZ') 
510  FORMAT(6F11.3) 

530  FORMAT(3F  11.  3) 
STOE 
2ND 


C - THIS  ROUTINE  CALCULATES  AND  PLOTS  THE  RZFRACTIVITY  AND 

C - REFRACTI VIT Y  GRADIENT  PROFILES 

C 
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IMPLICIT  E  EAL*8  (A-H,  0-  2) 

COMMON  /GNECCM/  CRH  ,  CRX  ,  CF.G  ,  DTR  ,  PH  AS  E,  I HETA  #C  ,  C  1 ,  C2 

1  ,  RKC  (51,30)  ,?  HI  (51, 3  C)  ,  PI ,  C  MF  ,  C  N  A’J  T,  DR  NCI 

2  ,  FRQ,TPB,  ABSR  H  ,  NHV,  NSL,  NF.M  S,  JAR2A,  NDATA 
COMMON  /T NO COM/  RA  ,  D ELX  ,  A ES ,  AHO  , SLO  (5  1)  ,  S Xfl AX,  X FINAL 

1  ,  J  R  A  Y 

CCMMCiJ  /THE  COM/  VX  (30)  ,  YfilN,  YMAX,  XEIY,  YDI7 

1  ,H  (51 ,3  0)  ,G(51 , 30)  ,  HN  (100)  , RN  ( 1 00)  , S  (5 1 , 3 0) 

2  ,  JPIT,  KEEF  ,XG R  A  E, K? 1CT, JPL  CT , IPLACE, JCASZ, NELO 

3  ,NREF,NEAY,NPRO,NPLCT, NGRAD, SCA1E 

COMMON  /FCRCCM/  DN  DZ  (1  00 )  ,R M  (10  0)  ,  DMDZ  (10  0)  ,  KMI  R 
1  , KKGFAD, KMIR, NMGRAE,HF (10  0) , EAICIL 

REAL  W,  Z 
INTEGER  I I ,LLL 

DIMENSION  X  (150)  ,Y  (1  50)  ,  E  (100) 

INITIALIZE  TEE  ARRAYS  X  AND  Y 
DO  5  1=1,100 
X (I) =0. 0 
Y  (I) =0. 0 

5  CONTINUE 

PRINT  HEALING  IE  THERE  IS  A  PRINTOUT 
FOP.  MODIFIED  INDEX  OF  R  SFEACTI VIT  Y 
IF  (JPLOT-5)  10,6,8 

6  CONTINUE 

IF  (KMIR)  18,18,7 

7  PRINT  130 
GO  TC  18 

8  CONTINUE 

IF  (KMGRAD)  18,18,9 

9  PRINT  1 3  5 
GO  TC  18 


C  PRINT  HEAPING  IF  THERE  IS  A  PRINTOUT 

C  FOR  EEFRACTIVITY 

10  CONTINUE 

IF  (JPLOT-1)  11,11,14 

11  CONTINUE 

IF  (KREF)  18,18,12 

12  PRINT  1 1 0 
GO  TO  18 

14  CONTINUE 

IF  (KGRAD)  18,18,16 
16  PRINT  120 
C 

C  SET  UP  INITIAL  CONDITIONS  TO  CALCULATE  REFR ACTIVI 

C  AND  REFRACTIVITY  GRADIENT  VERSUS  ALTITUDE 

18  CONTINUE 

IF  (NPRO-  1)  20,20,  30 
20  I D=  5  1 

DEL  H= (YMAX-YHIN) /5  0. 

A  =  Y  KAX+DELH 
GO  TO  40 
30  ID=  NEATA 
L  =  0 

a=o 

IF  (HN(I)-YMAX)  33  ,  5  8,  58 
33  CONTINUE 

IF  (HN(I)-YHIN)  20,3  5,  3  5 
35  DELH=(YHAX“HN  (1))/25. 

A=Y  MAX+  DELH 
ID=  2  5 
M  =  I E  -  1 
C 

C  SET  UP  A  LOOP  TO  PRINT  AND  PLOT  THE  FREE  SPACE  AN 

C  EXPONENTIAL  PROFILES . 
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40  CONTINUE 


DO  41  1=1,13 
A=A-EELH 

X  (I)  =C1  *3EXP  (-C2*A) 

41  Y (I) = A 

A  = Y  N AX+  D  EIK 
DO  55  1=1, ID 
A=A-ESLH 

IF  (JPLGT-1)  42,42,400 
400  CONTINUE 

IF  (JPLCT-5)  46  ,402,41  0 

PRINT  THE  EEFRACTIVIIY  PFCFILE 

42  X  1  =  X  (I) 

'/  1  =  Y  (I) 

IF  (HREF)  44,44,43 

43  PRINT  140, 1,  Yl, XI 

44  CONTINUE 

IF  (NREF )  54,54,50 

PRINT  THE  REFR ACTIVITY  GRADIENT  PROFILE 
46  CONTINUE 

IF  (I-ID)  47,54,54 

4  7  X  1  =  3281. 0*  (X  (I)  -X(I  +  1)  )  /  (Y  (I)  -Y  (1+  1)  ) 

Y  1  =  (Y  (I)  +  Y  (1  +  1)  )/2. 0 
IF  (KGRAI)  45,45,43 

48  PRINT  140,1, Y1, XI 

49  CONTINUE 

IF  (NGRAC)  54,  54,5  0 

402  GO  TC  54 
410  UC  TC  54 


c 


PLOT  THE  SET  F.  ACTIVITY  ANT  RE  ?F  A  CTI  VIT 1'  GP.  AD  IE 
C  ?  PC  f  ILE  5  . 

50  X1=X1/XDIV 

Y  1=  (Y  (I)  -  YMIN)  /  YDI  V 

IF  (1-1)  £1,51,52 

51  CALI  PLOT  (5NGL  (XI)  ,  SNGL  (Y1)  ,3) 

52  CALI  PLOT  (SNGL  (XI)  ,  SNGI  (Yl)  ,2) 

IF  (JPLOT-5)  413,54,53 

4  13  CONTINUE 

IF  (JPLGT-1)  54,54,53 

53  Y2=  (Y  (1+  1) -YKIN)/YDI7 

CALI  PLOT  (SSGI  (XI)  ,  SNGL  (Y2)  ,2) 

54  CONTINUE 

55  CONTINUE 

IF  (NPRO-1)  60,  80,  56 
C 

C  SET  OP  A  LCCP  TC  PRINT  AND  PLOT  THE  PISCE-WIS 

C  LINEAR  PROFILES. 

56  L=IE 
58  L?=0 

DO  75  1=1, NEAT A 
IF  ( H  N  (I)  -  YU  AX)  60,60,  74 
60  CONTINUE 

IF  (HN(I)-YNIN)  74,61,61 
6  1  L  =L  +  1 
L?=IE+1 
:i  =  M  + 1 

X  (L) =RN  (I) 

Y  (L)  =HN  (I) 

X  (1  +  1)  =EN  (1+  1) 

Y  (L  +  1)  =  H N  (1  +  1) 

Z  =  HF ( I )  /  1000.0 
E  (L)  =  X  (I)  +  (  157.  0  *  Z) 
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C  NOTE  THAT  Z  (1)  =  EM  (I)  IE  DATCTI  =  0, 

C  AND  E  (L)  HAS  NC  EE  1  ATI  CN  TO  EM  ( I)  IF  DATCTL  =  1 

IF  (JPLOT-1)  62,62,65 
C 

C  PRINT  THE  REFRACTIVITY  PROFILE 

o2  X  1  =  X  (L) 

Y  1  =  Y  (L) 

IF  (KRSF)  64,64,63 

63  PRINT  140,L,Y1,X1 

64  CONTINUE 

IF  (NREF)  74,74,70 
C 

C  PRINT  THE  REFRACTIVITY  GRADIENT  PROFILE 

65  CONTINUE 

IF  (JPLOT-2)  1  00,  1  00,200 

100  CONTINUE 

IF  (CATCTI)  66,66,  101 

101  XI  =  DNDZ  (I) 

Y  1  =  Y  (L) 

IF  (KGRAC)  65,69,68 

66  CONTINUE 

IF  (I-NDATA)  67,74,74 

67  X  1  =  3  281 . 0*  (X  (M)-X  (H+  1)  )/  (Y  (M)-Y  (M+  1)  ) 

Y  1  =  Y  (  L) 

IF  (KGRAE)  65,65,68 
63  PRINT  140,M,Y1,X1 
69  CONTINUE 

IF  (NGRAE)  74,74,70 
C 

C  PRINT  THE  NOTIFIED  INDEX  OF  REFRACTIVITY  PROFIT 

200  CONTINUE 

IF  (JPLOT-5)  201  ,20  1  ,  300 

201  CONTINUE 

IF  (CATCTI)  204  ,204,202 
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202  X  1  =  EM  (I) 
i  1  =  Y(L) 

IF  (KMIR)  210,210,208 
204  CONTINUE 

IF  (I-N  DATA)  206,74,74 
206  XI  =  E  { L) 

Y  1  =  Y  ( L ) 

IF  (KMIR)  210,210,208 
2C8  PRINT  140, I, XI, XI 
210  CONTINUE 

IF  (NMIR)  74,74  ,70 

C  PRINT  THE  MODIFIED  INDEX  CF  REF  RACTIVIT  Y  GRADIENT 

C  PROFILES. 

300  CONTINUE 

IF  (EATCTI)  204,304,  302 
302  XI  =  DMDZ  (I) 

Y 1  =  Y  (L) 

IF  (KGRAD)  312,312,310 
304  CONTINUE 

IF  (I--NDATA)  306,74,74 

306  LL  =  I  +  1 

IF  <IL  .IE.  NDATA)  GO  TO  308 

307  CONTINUE 

DMDZ  (I)  =  (E  (L)  /  HF  (I)  )  *  1000.0 
GO  TO  309 

308  W  =  RN(LI)  +  (157.0  *  ( HF  (LL)  /  1000.0)  ) 

DMDZ  (I)  =  (E  (L)  -  W)  /  ( E  F  (I)  -  HF  (LL)  )  *  1000.0 
3C9  X  1  =  DMDZ  (I) 

Y  1  =  Y  (  L) 

IF  (KGRAD)  212,212,310 
310  PRINT  140,1,11 , XI 
312  CONTINUE 

IF  (NMGRAD)  74,74,  70 
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c 


C  PLOT  THE  BEFEACTIVITY,  R  EFRACTI  VIT  Y  GRADIENT, 

C  MODIFIED  INDEX  CF  EE  FE  ACTIVITY,  AND  MODIFIED  INDEX  OF 

C  EEFEACTIV  IT  Y  GE  ADI  ENT  PSCFILES. 

70  X 1=  X 1/X D I V 

Y  1  =  (Y  ( M)  -  YM I  N)  /  YDI V 
IF  (IP-  1)  71 ,7 1 ,72 

71  CALI  PLOT  (SNGI  (XI)  ,SNGI  (Y1)  ,3) 

72  CALL  PLOT  (SNGI  (XI  )  ,  SNG  I  (Y1)  ,2) 

IF  (J PLOT-1)  74,74, 73 

73  Y2=  (Y(M+1)-YPIN)/YDIV 

CALI  PLOT  (SNGI  (XI)  ,  SNGI  (Y2)  ,2) 

74  CONTINUE 

75  C  CNTIN'JE 
C 

C  POSITION  THE  PEN  IF  THERE  HAS  BEEN  A  PLOT 

80  CONTINUE 

IF  (J PLOT -5)  88,82, 34 

82  CONTINUE 

IF  (NMIE)  56,96,94 
84  CONTINUE 

IF  (NMGRAD)  96,56,  94 

83  CONTINUE 

IF  (JPLOT-1)  90,90,92 
90  CONTINUE 

IF  (NPEF)  96,56,94 
92  CONTINUE 

IF  (NGRAD)  56,  96,94 
94  CALI  PLOT  <0- ,0- ,-999) 

CALI  PLOT  (2. 0,2. 0,-3) 

CALI  FAC  T  C  E  (S  NGL  { S  C  AL  E)  ) 


96  CONTINUE 
106  RETURN 


110  FORMAT  (///23X, 'REFRACTI  VITY  PR  CFIIE  '/I  OX  ,  •  I  *  ,  3X, 

1  *  ALTITUDE  ' , 3X,  ’  RE?  F.  ACTI  V ITY  • /2 1  X,  •  (FT)  »  ,  1  1 X , 

2*  (N-UNITS) • ) 

120  FORMAT  (///  1  SX ,  '  RE  FR  AC  TI  VII Y  GRADIENT  PRO  FI LE '  /I  0  X,  .. 
1*  I'  ,8X,  ’  AITITUDE',7X,'R£FF  GRAD  IENI  '/21  X,  *  (FT)  *  ,  9  X,  .. 

2  *  (N- UNIT S/KM ) ' ) 

130  FORMAT  (/// 9  X,  '  MOD  I F  IE  D  INDEX  OF  REFRACTI  VITY  ', 

1PROFILE'/10X,'I',8X, 'ALTITUDE1 ,  5X , 'MOD  INDEX  OF  RSFR ' , 
2/21  X,'  (FT)  1  ,  11X,'  (M-UNITS)  ») 

135  FORMAT  {///I  OX, ’ MO DIFIE D  INDEX  CF  REFRACTI VITY  ', 

1  ’  GRADIENT  FRCF  HE’  ,/ 10  X,  » 1*  ,  3X,  *  ALTITUDE'  ,  7X  , 

2*  MOD  INDEX  OF  REFR  GRADIENT1, 

3/21  X,'  (FT)  '  ,9X,  ’  (M-D  NITS/KM)  •) 

140  FORMAT  ( 8  X,  1 3 , 2  (6X  ,  F  10  -  2 )  ) 

END 

C************************* ********************************** 

S  UBROUTINE  PKINTG 
C 

C - THIS  ROUTINE  SETS  UP  THE  INTEGRATION  OF  THE  F. AY 

C - TRAJECTORY  AND  TIME  OF  EROAGAT ION  EQUATIONS  FOR  EACH 

C - RAY.  WHEN  A  RAY  CROSSES  A  BOUNDARY  EETWEEN  THE 

C - THE  PIECE-WISE  LINEAR  SEGMENTS  OF  THE  REF  RACTI VIT  Y  .. 

C - PROFILE  CR  TEE  BOUNDARY  AT  THE  EARTH'S  SURFACE,  THE  .. 

C - RAY  EQUATIONS  ARE  INTEGRATED  TO  THE  BOUNDARY  BY  MEANS 

C - OF  A  VARIAELE  STEP  SIZE  INTERPOLATION  ALGORITHM.  THE 

C - RAY  EQUATIONS  ARE  THEN  RE-INITIALIZED  AT  THE  BOUNDARY 

C - AND  INTEGRATED  TO  THE  NEXT  BOUNDARY,  WHERE  THE 

C - INTERPOLATION  IS  REPEATED. 

C 

IMPLICIT  F  E AI*8  (A-H,  O-Z) 
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COMMON  /CNECCM/  CF.  H,CF.  X, CFG,  DTR, PHASE, THETA, C,C1,C2 

1  ,EHC  (5  1  ,30)  ,PHI  (51 , 3C)  ,  PI,  CMF,CNA'JT,  DEN OT 

2  ,FRv,TFfc, AESRH, NHV, NSL,  NR MS ,J AREA, HD AT A 
CCHHCN  /TKCCCH/  RA  ,  D  ELX  ,  AHS,  AHO  ,  ELC  ( 5  1)  ,EXMAX,XFINAL 

1  ,J  FAY 

CCHHCN  /THRCCM/  VX  ( 3  0)  ,  THIN,  YM  A  X,  XDI 7  ,  YDI  V 

1  ,H  (5  1  ,3C)  ,G  (5  1  ,  30)  ,  HN  {100)  ,  RN  ( 1  0  0)  ,  S  (5  1 , 3  0) 

2  , JPLT,KBEF,KGRAD, KP ICT, JPL CT , IPLACE, JCASE, NELO 

3  ,NREF, NRAY,NPRO,N PLOT, NGR AD, SCALE 

CCHHCN  /FCBCCH/  DNDZ  (1  00)  ,HH  (10  0)  ,  DMDZ  ( 10  0)  ,  KHIH 
1  ,  K  M  G  R  A  D  ,  N  M IR  ,  NHGRAE/HF  (10  0)  ,  DATCTL 

DIMENSION  TINT  (10)  ,  DELXX  (5)  ,  PINT  (10) 

DATA  DELXX,NEQ/100  00  .DO,  1  GOO .  D  0 , 1  0  0.  DO,  1 0  .  D  0 , 0 .  DO  ,  2/ 

C  SET  UP  INITIAL  CONDITIONS 

DO  5  1=1 , 10 
PINT  (I)  =0.0 
5  TINT  (I)  =0.0 
L  =  1 

ANG1E=0.0 
RHHAG=1.0 
S  TP  X  =  DEL  X 
VX  (I)  =C R X 
3  ( JRAY, L ) =CE  K 

S  ( JF  AY,  L)  =  I A  TA  N  (CR  G*  RA/  (  RA  +  CRH)  ) 

G  (JRAY,  L)  =  EI  NT  ( 1) 

F  HO  (JRAY, I)  =  F  HM  AG 
PHI  (JRAY,  I)  =  ANG IE 
YINT  (1)  =CFH 
YINT  (2)  =CRG 

V  A  LG  =  YI  NT  (2) 

V GT EHP=  PINT  ( 1) 

VHTEHP=  YINT  (1) 

V  XT  EM?=C  S X 
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X  =  C  R  X 
ICCDE=3 
C  11  =  0 

11=  1 

c 

C  CHECK  IF  A  RAY  TRACE  IS  TC  BE  MADE 

IF  (NRAY)  14,14,7 
C 

C  TEST  FOR  MAXIMUM  AND  MINIMUM  ALTITUDES 

7  IF  (CRH- Y MAX)  8,8,  12 

8  IF  (CF.H-  YMIN)  13,1  0,  10 
C 

C  ALTITUDE  EOUNDS  NOT  EXCEEDED.  POSITION  PEN  AT 

C  EMITTER  CCCRCINATES. 

10  Y  =  (CRH-YMIN)/YDTV 

CALI  PLOT  (SNGL  (X)  ,  S NG  L  (  Y)  ,  ICOD E) 

ICGDE=2 
GC  TC  14 
C 

C  MAXIMUM  AITITIDE  EXCEEDED.  POSITION  PEN  AI  UPPER 

C  LEFT  HAND  GF.AEH  CORNER. 

12  Y= (YMAX-YMIN) /YDIV 

CALI  PLOT  (SNGL  (X)  ,  S  NG  I  (Y)  ,  ICODE) 

GO  TC  14 
C 

C  MINIMUM  AITITIDE  EXCEEDED.  POSITION  PEN  AT  LOWER 

C  LEFT  HAND  CCFNEE. 

13  Y  =  0 . 0 

CALI  PLOT  (SNGL  (X)  ,SNGL(Y), ICODE) 

C 

C  FIND  WHICH  LAYER  THE  EMITTER  IS  IN 

14  IF  (GRAY  -  1)  15,15,22 

15  CONTINUE 


100 


DO  201=1 ,NCATA 
IF  (CEK-hS  (I) )  20/  20,16 
16  J  AR  E  A=I 
J  X  .1 1 R  =  I 
GO  IC  25 
20  CONTINUE 
GO  TC  25 
22  JAREA=JXHTE 
25  CONTINUE 

SET  UP  A  IOC E  FOR  INTEGRATION  Of  THE  ARRAYS  "YINT" 
AND  "PINT". 

DO  2  CO  1  =  2, IELACE 
11=11+1 


CALI  ROUTINE  TC  INTEGRATE  "YINT"  AND  "PINT" 

CALI  RK  (NEQ,CRX,STPX,  YINT, PINT) 

CHECK  WHICH  IAYER  THE  RAY  IS  IN 
100  HARE  A= J  AREA 

DO  120  J  =  1 ,  N I AT  A 
IF  (YINT  (l)-HN(J))  120,  1  20,  1  10 

110  K  AR  E  A=J 
GO  TC  125 
120  CONTINUE 

SET  UP  LAYER  IF  RAY  HAS  INTERSECTED  EARTH'S  SURFACE. 
IF  (YINT  (1) -AHS)  1  22,1  22,125 
122  RARE  A=N  D  AT  A  +  1 
125  CONTINUE 

CHECK  WHICH  IAYER  BOUNDARY,  IF  ANY,  HAS  BEEN  CROSSED 


FIRST. 


(JAREA-KAFEA)  140,145,130 


AN  UFPE  R  ECUNDAFY  HAS  3  E  EN  CROSSED 
130  3NDF.Y  =  HN  (JAEEA-  1) 

KAE  EA=JAR  EA- 1 
GO  1C  160 


A  LC  ViSR  BOUNDARY  HAS  3  EE  N  CROSS  ED 
140  3NDEY=HN  (JAEEA) 

RARE  A=J  A  EE A  +  1 
GC  TC  160 


NO  BOUNDARY  HAS  BEEN  CROSSED-  STORE  ARRAY  VALUES 
EVERY  TENTH  INTEGRATION  STEP. 

145  IF  (11-  10)  148,  146,  1  46 

146  11=0 
L  =  L  +  1 

VX  (1)  =C  F.  X 
H  ( J  R  AY,  L )  =Y  INT  (1) 

S  ( J R  A Y ,  L )  =D  ATA  N  (YINT  (2)  *  f  A/  (RA  +  YINT  (  1)  )  ) 

G (JFAY, L) =PINT  (1) 

RHO  ( JRA  Y  ,  I)  =RHMAG 
PHI  (JEAY  ,1)  =  ANG IE 
143  V ALG  =  YINT  (2) 

VXTEMP=C  F  X 
VHTE«P=YINT  ( 1) 

V  GTE  H.P=PI  NT  (1) 

5TP  X  =  DEL  X 


CHECK  IF  A  RAY  TRACE  IS  TC  BE  MADE 
IF  (NRAY)  200,200,  152 


A  FAY  TRACE  IS  TO  3E  MADE 
MINIMUM  AITITUDE. 


CHECK  FOR  MAXIMUM  AND 


152  IF  (YINT  (1)  — YMAX)  1  53,  152,150 

153  IF  (YINT  (  1) -Y.MIN)  156,156,154 


MAXIMUM  ALIITUDZ  NOT  EXCEEDED.  CALL  ROUTINE  TO  PLOT 
THE  FAY. 

154  X  =C  FX/X  D 1 7 

Y=  (YINT  (1)-YHIN)/YDIV 

CALI  PLOT  (SNGL(X)  ,  S  NG  L  (  Y)  ,  ICQD  E) 

I  CO  EE=2 
GO  TO  158 

MAXIMUM  OF  MINIMUM  ALTITUIE  EXCEEDED.  TURN  OFF  PLCTTE 
156  IC0EE=3 

158  IF  (JARE A-K AREA)  170,200,170 

THE  LAYER  ECUNEARY  HAS  3EEN  FOUND.  SET  UP  FOR  LINEA 
INTERPOLATION  SCHEME. 

16  0  D  X  TC 1  =  0 . 

H?EME=YINT  (1) 

C  R  X  =  V  XT  E  M  E 
Y  INT  (1)  =  VRTEMP 
YINT  (2)  =  V ALG 
PINT  (1)  =  V GTE MP 


SET  UP  VAHIAELE  INTEGRATION 
TO  THE  3C ONE  ARY. 


STEP  SIZE  AND 


INTERPOLA 


TE 


DO  165  I J  K=  1 , 5 
X  X  =  -  Y I N T  ( 2)  /  YI  NT  (3) 

X  YCHK=XX*XX-2.  *  (YINT  (1)  -ENDEY)  /YINT  (3) 
IF  ( XYC HK.L1.0.)  X YC HK=  0  . 

Y  Y=  ES^RT  (XYCRK) 

CHGX  =  XX  +YY 


IF  (XX.GT.YY)  C  HGX=  X  X-  Y  Y 

IF  (IJK.GT.4)  CHGX=  (3NDR  Y-YINT  (  1)  ) /YINT  (2) 


5D=ST?X-DXTCT 

IF  (  (CHGX.LI.O.)  .0?,.  (CHGX.GE.SD  J)  CHGX  = 

1  ( 3  T  E  X-DXTCT)  *  ( 3ND  E  Y-YI NT  (  1)  ) / ( HTEMP- YI NT  ( 1 )  ) 

CHGX=CHGX-DZIXX  (IJK) 

IF  (CHGX.IE.O.)  GO  TC  165 
C 

C  CALI  ROUTINE  TC  INTEGRATE  TO  THE  3  CU  N  DARE 

CALI  RK  ( N EQ,  C F. X, CHGX,  YI NT , ? IIIT ) 

DXTCT=DXTCT+CHGX 
165  CONTINUE 
C 

C  CHECK  IF  A  RAY  TRACE  IS  TC  BE  HADE 

IF  (NRAY-1)  170,152,  152 

r 

C  CHECK  IF  FAY  HAS  INTERSECTED  EARTH'S  SURFACE 

170  IF  (KAREA- (NEATA+1)  )  190,180,  18  0 

C 

C  RAY  HAS  CEC5SED  ZERO  ALTITUDE  BOUNDARY.  FIND  INC  IDEM 

C  GRAZING  ANGLE. 

180  T  HET  A=D  A  B  S  (E  AT  A  N  (YI  NT  (2 )  )) 

C 

C  CALL  ROUTINE  TC  CALCULATE  CCHPLEX  SCATTERING 

C  COEFFICIENT. 

CALL  SCATT 

C 

C  SET  UP  CCHELEX  SCATTERING  COEFFICIENT 

RHMAG=RHHAG*AESRH 
A NGLE=ANGLE+PHASE 
C 

C  ADD  KULTIFATK  RAY 

DO  185  L L  =  2 ,  10 
185  YINT  (LL)  =-YINT  (LL) 

KAREA=MDATA 

C 

104 


SET  UP  TC  INTEGRATE  FROM  THE  BOUNDARY  TO  THE 
"OELX". 

Isu  JARZA=KAFEA 
VXTEM?=CFX 

V  HTE  PP=  Y I  NT  (1) 

V  ALG  =  YI  N  T  (2) 

VGTER?=?INT  (1) 

3TPX=STPX-DXTCI 

CALL  ROUTINE  TC  INTEGRATE  FF.O.T  THE  BOUNDARY 
NEXT  "DELX". 

CALL  FK  ( NE ,C R  X,  5 T?  X,  YI NT, PINT) 

CHECK  FCE  MOPE  EOUNDARY  CROSSINGS 
GO  TC  100 
200  CONTINUE 
ICODZ=3 

IF  ( NR A  Y )  230,230,  2  1  ^ 

210  CALL  PICT  (SNGL  (X)  ,S  NGL  (Y)  ,  ICODE) 

CHECK  IE  THIS  IS  A  PLOT  CF  THE  LAST  RAY.  IF 
POSITION  THE  PEN  FOR  THE  NEXT  PICT. 

IF  (JRAY-NELC)  230,2  20,  2  20 
220  CALI  PLOT  (0.,0  .,-999) 

CALL  PLOT  (2. 0,2. 0,-3) 

CALL  FACTCF  (SNGL  (SCALE)  ) 


230  CONTINUE 


S  LIE  EC  UT I N  E  R  K  ( N  ,XM,  H,  Y,  r) 

C 

C - THIS  ROUTINE  INTEGRATES  THE  FA  x  TRAJECTORY  AMD  TIME  OF 

C - PROP  AGATICN  DIFFERENTIAL  EQUATIONS  3Y  MEANS  OF  A 

c - FOURTH  CFIE5  RUNGS-KUITA  ALGOEI IHM  USING  AN 

c - -INTEGRATION  STEP  SIZE  CF  H  =  XDELTA/10,  WHERE  XDELTA  . 

0 - IS  TEE  DISTANCE  INTERVAL  ALONG  THE  EARTH'S  SURFACE  .. 

C - FOR  FEINTING  ANE  PLOTTING  THE  HEIGHT  GAIN  CURVES,  OR 

C - H  =  CHGX  REZEE  CHGX  IS  A  VARIABLE  STEP  SIZE  SET  BY  THE 

C - INTERPOLATION  ALGORITHM  IN  SUBROUTINE  RKINIG. 

C 

IMPLICIT  EE AI*9  (A-H,  C-Z) 

DIMENSION  Y  (10)  ,P(10)  ,YDCT  (10)  ,  EDCT  (10)  ,Q  (11,4) 

1,  R  (  10,4)  ,  YN  (10)  ,PN  (1  0) 

C 

C  SET  UP  INITIAL  CONDITIONS 

DO  5  1= 1 , N 
YN(I)=Y(I) 

5  ?N  (I)  =P  (I) 

C 

C  SET  CP  A  LOOP  TO  INTEGRATE  THE  DIFFERENTIAL  EQUATIONS 

DO  6C  L= 1  ,4 

c  .  . 

C  CALL  ROUTINE  TO  SET  UP  TRr  DIFFERENTIAL  EQUATIONS 

CALI  ATMOS  (  Y,  P  ,  YDCT  ,?  DOT) 

C 

C  INTEGRATE  TEE  DIFFERENTIAL  EQUATIONS 

GO  TO  (10,2  1,30  ,40  )  ,  L 
10  DO  IF  K  =  1  #  N 


Q  ( K  ,  I )  =  H  *  Y  P  C  T  (  K ) 

15  Y  (K)  =y:j  {F)*Q  (R,I)/2 


a  ( 1 , 1)  =  H  *  FD  Cl  ( 1 ) 

P  (1)  =PN  (  1)  +  R 
X=XN+H/2 . 

GO  1C  50 
20  DO  25  K  =  1  ,  N 

Q  (K,L)  =H*YDC1  ( K) 

25  Y  (K)  =  YN  (K)  +5  (K,I)/2. 

H  (1  ,1)  =H*EDCT  (1 ) 

?  (1) =PN  (1)+R  ( 1  ,  L)  /2  . 

X  =  X  N  +  H/2 . 

GO  1C  50 
30  DO  35  K= 1 ,N 

Q  (K  ,  1)  =  H*  YDC 1  (K) 

35  Y  (K)  =  YN  (K)  +  Q  (K  ,1) 

R  ( 1  ,  L)  =  H*  EDCT  (1) 

P  (1)=PN(1)+R  (1,1) 

X  =  X  N  +  H 
GO  1C  50 
40  DO  45  K  =  1  ,  N 

2  (K  ,1)  =  H *  YDC1  (K) 

4  5  Y  (K)  =  YN  (K)  +  (Q  (K,  1)  +2.0*Q  (K,  2)  +2 . 0*Q  (K,3)  +  Q  (K,  4)  )  /6.  0 
S  (1  ,1)  =H*FDCT  (1) 

P  (  1)  =PN  (  1)  +  (E  (1 ,  1)  +2 .0*R  (1, 2)  +2 .  0*E  (1,3)  +  R  (1 ,4)  )  /6.0 
X  N=  X  N  +  H 
50  CCN1INUE 
60  C  CN IIN'JE 

CALI  ROUTINE  1C  FIND  THE  VALUES  OF  THE  DIFFERENTIAL 
EQUATIONS  AT  THE  END  OF  THE  INTEGRATION  STEP. 

CALL  AT  .'1C  S  (  Y,  P  ,  YDOT  ,P  DOT ) 


STORE  THE  NEW  DERIVATIVE  CF  THE  RAY  SLOPE 
Y  (3)  =  Y D C T  (2) 

RETURN 


o: 


*******3**********  ***** 


£*********#**$**$**************#******************$********* 
S  UBROUTI  NE  ATMOS  (  Y,  P,  YDCT,  PDOT) 

c  •  « 

c - THIS  HOOTINE  computes  atmospheric  refractivity  at 

C - THE  RAY  ATTITUDE  AND  SETS  UP  THE  RAY  DIFFERENTIAL 

C - EQUATIONS  . 

C 

IMPLICIT  REAL*8  (A-H,0-Z) 

CCMMCN  /CNECCM/  CRH  ,  CRX  ,  CEG,  DTR  ,  PHAS  E,  THE  7A  ,  C,  Cl  ,  C2  . 

1  ,RHO  (51  ,30)  ,PHI  (51 , 30)  ,PI,  CMF,CNA'JI,  DSNOT 

2  ,FRQ,TPK, AESRH, NHV, NSL,  NE M  S  ,  J  AREA,  ND  ATA 

CCMMCN  /TKOCCM/  RA  ,  DELX  ,AES,AHO,ELC(51)  ,  EXMAX,  XFINAL  . 
1  /JR AY 

CCMMCN  /THBCCM/  VX  (3  0)  ,  Y  MIN,  YMA  X  ,  XDIV  ,  YDI V 

1  ,H  (5  1,30)  ,G(51, 30)  ,  HN  (100)  ,RN  (100)  ,S  (51 ,30) 

2  /JPIT/KREF/KGRAE,  KPICT,  JPLCT/IPLACE,  JCASE,NELO  . 

3  ,NREF,NEAY,NPRO,NPLCT,NGRAC, SCALE 

COMMON  /FCRCCM/  DNDZ  (1  00)  ,RM  (10  0)  ,  DHDZ  ( 10  0)  ,  KMIR, 

1  ,  KMGRAD,  KMIR,  NMGRAE/HF  (10  0)  ,  DATCIL 

DIMENSION  Y  (10)  ,P(10),YDOT(10)  ,  PDOT  ( 1  0) 

C 

C  SET  UP  INITIAL  CONDITIONS 

C  H=  Y  ( 1 ) 

CG=  Y  (2) 

C 

C  TEST  FOR  THE  APPROPRIATE  ATMOSPHERIC  MODEL 

IF  (NPRO-1)  10,10,20 
C 

C  FREE  SPACE  AND  EXPONENTIAL  MODELS 

10  R  EFB  =  C1  *  E  EXP  (-C2*CH) 

DLNEH=-C2*EEFR*1.0E-06/ (REFR  +  1.  0E-06  +  1.0) 


GO  TC  IOC 


THE  ATMO  SPHERE  IS  STRATIFIED.  SELECT  THE  APPROPRIATE 
MODEL. 

20  IF  (JA REA-1)  50  ,50,6  0 


ALTITUDE  IS  ABOVE  THE  HIGHEST  R  EFR  ACTIV IT  Y  PROFILE  . 
DATA  POINT.  USE  AN  EXPONENTIAL  MODEL  WHICH  FITS  THE 
DATA. 

50  SEFR=C1*E  EXP  (-C2+CH) 

0LNEH=-C2*EZFR*1.0E-06/(BEFR*1.  0E-06  +  1.0) 

GO  TC  100 

ALTITUDE  IS  EZLCW  THE  HIGHEST  R  EFR  ACTIV IT  Y  PROFILE  . 
DATA  POINT.  USE  A  PIECE-WISE  LINEAR  MODEL. 

60  SLC  FE=  (F.N  (JAREA-1)  -R  N  ( J  AREA)  )  /  (  HN  (JAREA-1)  -HN  (JAREA)  ) 
3  =  RN  (JAF.EA) -SLOPE*HN  (JAREA) 

• 

COMPUTE  R  EFRACI IVITY  FOR  THE  PIECE-WISE  LINEAR  MODEL 
30  R  EFR=SLO  P  E*CH  ♦  B 

DLNCK= (SLCFE*1 . 0E-06)/ (REFR  +  1.  0 E-06  +  1.0) 

• 

COMPUTE  THE  DERIVATIVES  FOR  RAY  TRACES 
00  E  AD  =  EA*C  H 

YDOT  (2)  =ELNDH*  (  (EAD/RA)  **2  +  CG**2)  +  (2.  0/RAD)  *CG**2 
1  +  RAD/R A  **2 
Y  CO  T  ( 1 )  =  C  G 

• 

COMPUTE  THE  DERIVATIVE  FCE  TIME  OF  PROPAGATION 
PDOT  (1)  =RAE*  (REFR*1. 0E-06  +  1.0) 

1*DS£ET{1 .0+  (FA*CG/R  A  D)  **2)/(C*RA) 

RETURN 


SUBROUTINE  SCAT! 

C 

C - THIS  ROUTINE  CALCULATES  TEE  COMPLEX  SPECULAR 

C - SCATTERING  COEFFICIENT  FOR  WAVEFRONT  REFLECTION  FROM 

C - SMOOTH  OR  ROUGH  LAND  AND  SEA  SURFACES. 

C 

IMPLICIT  F£AI*8  (A-H,  O-Z) 

COMMON  /CNECCM/  CRH  ,  OR  X  ,  CRG,  DTR  ,  P  HAS  E,  I  HET  A  ,  C ,  C  1  ,  C2 

1  ,  RHC  (5  1  ,30)  ,PHI  (51  ,  3G)  ,PI,  CHF,CNAUT,  D3NCT 

2  ,FRQ,TFR, ABSRH,  NH V , NSL, NR M S , J AREA, ND  AT A 
DIMENSION  EPS  (3,2)  ,  SIGMA  (3,2)  ,  D  ELHSL  (9,  2) 

REAL  NI,  NR 

DATA  EPS,  IC/ 80.  DO,  6  9 .DO, 6 5. DO, 4 .D0,1 0.D0, 30. DO, 0/ 
DATA  SIGMA/4. 3D0,6. 5  CO, 1.6D1,1.  D-4, 1 . 6D-3 , 1 . 0D-2/ 
DATA  DELHSL/0.  DC,  2. 0D-  1, 6.0D-1,  1 .  1  DO  ,  1.  7D  0 , 2 . 6D  0 
1,  4.  3D0, 8. 6D  0,1. 29D  1 , 0.  D0#  9.  DO,  3  0.  DO,  5  6.  DO,  11  2.  DO 
2, 2 14. DO, 4 2 9. DO, 128 8. DO 
3,  2.  1 46D 3/ 

C 

C  SET  UP  INITIAL  CONDITIONS 

IF  (IC)  10,10,50 
10  IC=1 

J  J  =  NRMS  +  1 

IF  (NSL-1)  12,12,30 
12  K  K=  1 

IF  (FRQ-  15000C0000.)  14,  16,  16 

14  11=1 

3  0  TO  4  0 

16  IF  (FRQ-5OC0COOOOO.)  18,20,20 
18  11=2 

GO  TO  40 


20  11=2 

GO  1C  40 
30  KK= 2 

1 1=  NSL-  1 
C 

C  CALCULATE  THE  COMPLEX  DIELECTRIC  CONSTANT 

40  3fi=EFS  (II, KK) 

EI  =  60. 0*  SIGMA  (II,  KK)  *C/  (CHF*FRQ) 

A  =  D  £  CRT (EE**2+EI**2) 

A LP K A=D A T AN2  (EI,ER) 

DELTA H=DEIHSL  ( J  J  ,  K  K) 

C 

C  CHECK  WHICH  POLARIZATION  IS  BEING  USED 

50  IF  (NHV-1)  60,60,70 

C 

C  HORIZONTAL  POLARIZATION 

60  N  R=  C  SIN  (TRETA)  -DSQRT  (A)  *DCOS  (ALPHA/2.) 

NI=ESCRT  (A)  *ISIN  (ALPHA/2.) 

D  R=  E  SIN  (THETA)  ♦DSQRT  (A)  *ECCS  (ALPHA/2.) 

DI  =  ESQRT  (A)  * E S I N  (ALPHA/2  .) 

G  R=  (KR*DE-NI*DI)/(DB**2*EI**2) 

G  1=  (NR*DI*NI*DR)/(DR**2  +  EI**2) 

GO  TO  30 
C 

C  VERTICAL  POLARIZATION 

7  0  NR  =  E  SQRT  (A)  *ECCS  (ALP  HA/2  .)  *DSIN  (THETA)  -  1.  0 
NI=  DSQRT  (A)  *  ESIN  (AL  P  HA/2  .  )  *  DSTN  (THETA) 

DR=ESQRT  (A)  +000  5  (ALPHA/2  .)  *DSIN  (THETA)  ♦  1.  0 
DI=D£QRT  (A)  *ESIN  (AL P HA/2  . )  *DSIN  (THETA) 

GR=  (NF.*DR+NI*DI)/(DR**2  +  Ei**2) 

31=  (NE*DI-NI*DF.)/(DR  **2+DI**2) 

C 

C  CALCULATE  THE  COMPLEX  FRESNEL  REFLECTION  COEFFICIENT 

30  A3SRH=DSQ?T  (GR**2+GI**2) 


1  1  1 


P  H  A  S  H=  DA  T  AN  2  (-GI,3F.) 


C 

C  CALCULATE  TEE  TCTAL  SPECULAR  REFLECTION  COEFFICIENT  .. 

A  BSF,  H=A  BSFH*EZXE  (-0-  5*  (4  -  0*PI*D  ELT  AH+DSIN  (THETA) 

1  *FRC/C)  **2) 

RETURN 

END 

C 

£$*$***  **********-****$**$$*S***$*:fr4< *  *  *****$*$*****$♦£****$£# 

SUB  ROUTINE  HGAIN 
C 

C - THIS  ROUTINE  COMPUTES  FIELD  STRENGTH  AND  POWER  DENSITY 

C - HEIGHT-GAIN  FUNCTIONS  (ONE-WAY  FIELD  STRENGTH  AND 

C - POWER  DENSITY  NORMALIZED  TO  FREE  SPACE  VALUES)  AT 

c - "X DELTA"  DISTANCE  INTERVALS  ALONG  THE  SURFACE  OF  THE 

C - EARTH.  HEIGHT-GAINS  ARE  CETAINED  FROM  CALCULATIONS  AT 

C - 200  "WINDOWS",  OR  ALTITUDE  INCREMENTS,  EXTENDING 

C - VERTICALLY  EETWESN  THE  HIGHEST  AND  LOWEST  RAYS  WITHIN 

C - THE  SPECIFIED  PIOT  LIMITS. 

C 

IMPLICIT  FEAI*8  (A-H,0-Z) 

COMMON  /CNFCCM/  CRH , CRX , CRG, DTR , PHAS E, T HE TA , C, C 1 , C2  .. 

1  ,RHC  (51,30), PHI  (51, 3C)  ,PI,  CMF,CNAUT,  DENOT 

2  ,FRC,IPN, AESRH, NHV, NSL, NRM S , J A EE A, NDATA 

COMMON  /TWOCCM/  BA  ,  DELX  ,  AHS,  AHO  ,ELC  (5  1)  ,  EXMAX,  XFIHAL  . 
1  ,  JR  A  Y 

COMMON  /THRCCM/  VX  ( 3  0)  ,  Y  MIN,  YM A  X  ,  XDI V  ,  YDI  V 

1  ,H  (51, 3C)  ,G(51,30)  ,HN(100)  ,RN  (100)  ,S  (51,30) 

2  ,  JPLT,  KREF  ,KGRAD,  KPICT,  JPLCT,  IPLACE,  JCASE,  NELO  . 

3  ,NREF,NEAY,NPRO,NPLCT,NGE AD, SCALE 

COMMON  /FCECCH/  DNDZ  {1  00)  ,RM  (10  0)  ,  DMDZ  (  10  0)  ,  KMIR 
1  ,  KMGBAD,  NMIR,  NMGRAD,KF  (10  0)  ,  DATCTL 


r 

K 

m 

>■ 

_  • 

a 
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DIMENSION  L  (51)  , TD  113(30), ESIG(  30)  ,PSC(30) 
1,  TSIG  (3  0)  ,ESIGT  (30)  ,  PSCT  (30)  ,TSIGT  (3  0) 

DATA  ZDIV/40.0D0/ 


C 


C 

c 

c 


SET  UP  INITIAL  CONDITIONS 

DO  5  1=1/30 

TDIFF(I)  =0-0 

ESIG  (I)  =0.0 

PSC  (I)  =0 . 0 

TSIG  (I)  =0.0 

ESIGT  (I)  =0.0 

PSCT  (I)  =0.0 

TSIGT  (I )  =0. 0 

5  CCNTINOE 

DO  6  1=1,51 
L (I) =0. 0 

6  CONTINUE 

3  1  =  E A  +  A  HC 
XN=E£LX/XDIV 

Y  N=  YMAX/ YEIV 

Y  S=AHS/YDIV 
N  P  =  2  0  0 


SET  UP  A  LOOP  TO  CALCULATE  RELATIVE  FIELD  STRENGTH  OR 
POWER  DENSITY  AT  EACH  INCREMENT  OF  DISTANCE. 

DO  1000  K  =  2 , IPL  ACE 
H LO  =  YMA  X 


U=VX  (K)  /RA 
C 

C  FINE  HIGHEST  AND  LOWEST  RAYS 

DO  ICO  1=1,  NELC 
H  LO  =  EM  IN  1  (HLC/H  (I,  K)  ) 

100  CONTINUE 

H  H  I  =  D.MI  N  1  (H  ( 1,  K)  ,  Y  M  A  X) 


1  1  3 


HDIFF=HHI-HLC 

HTB  =  DA3  5  (  F  ( 1  , K)  -K  ( 2 ,  K)  ) 


EXIT  IF  AIL  SAYS  EQUAL  OS  EXCEED  THE  MAXIMUM  ALTITUDE 
IF  (HLO-YMAX)  150,  1  100  ,  1  100 

SET  U?  WINDOW  SIZE  AND  WINDOW  POSITIONS  TO  EXCLUDE  THI 
HIGHEST  AND  LOWEST  RAYS. 

150  Vi  INECW=HDIFF/1S0. 

H NO K=HHI- WINDOW/ 10. 

H  FI N  AL=H1 C  +  W 1NDCW/ 10  . 

RWP=HNOW-E FINAL 


FINE  THE  ALTITUDE  INCREMENT  FOR  EACH  NEW  WINDOW 
POSITION. 

DH=RWP/FICAT  (NP) 

HNOW=HNOW+DH 

SET  UP  A  LOOP  TO  POSITION  THE  WINDOW  EVERY  ”DH"  FEET 
IN  AITITUEE. 

DO  9C0  J=  1 , N  F 
H  N  0  W  =  HN  0  W  -  D  H 


CALCULATE  THE  UPPER  ANGLE  LIMIT  SUETENDED  BY  THE 
WINDOW  IN  FREE  SPACE. 

R  2=  R A  +  HNC  W 


D 5Q  =  R 1  *  *2  +  R2 *  *2  -  2. 0 *R  1*R2*DCOS 
E  1  =  EARCOS  (  (DSQ  ♦  {R  1  *DTA  N  (U)  )  **  2  - 
1/  (2.0*R  1  * E T A N  (U)  *DSQRT  (DSQ  )  ) 

IF  (R1/DCCS  (U)  .  GT.  R  2)  E1=-E1 


(UJ 


(R  1/DC  OS  (rJ)  -R2)  **2  ) 


CALCULATE  THE  LCWER  ANGLE  LIMIT  SUBTENDED  BY  THE 
WINDOW  IN  FREE  SPACE. 

R2=EA+H NOW- WINDOW 


DSQ  =  E1**2  +  52**2  -  2.  0  *F  1*R2*DC0S  (U) 

E2=  E ARCCS  (  (DSQ  +  (  5  1  *D  TA  N  (U )  )  **  2  -  (R  1/ECOS  (rJ)  -F.2)  **2  ) 
1/  (2-0*R  1  *ETAN  (U)  *DSQ5T  (DSQ  )  ) 

IF  (R1/DCCS  (U)  .  GT.  R  2)  E2=-E2 
SDIFF=DSQFI  (  DA  E  S  (E1-E2)  /  HR) 

CALCULATE  THE  ELEVATION  ANGLE  OF  THE  EMITTER  AT  THE  . 
WINDOW. 

R  1  =  R  A+  HNC  W-VI INDGW/2 . 

R  2=RA+A  HC 

D  SQ  =  R 1  *  *  2  +  F.  2  *  *  2  -  2.  0*F  1*R2*DCOS  (U) 

AQ  =  EARCOS  (  (ESC  ♦  (R  1  *DTA  N  (U)  )  **  2  -  (F.  1/DCOS  ( TJ)  -R2  )  **2 

1/  (2 . 0*R  1*EIAN(U)*D5QRT  (DSQ  )  ) 

IF  (S1/DCCS  (0)  .GT.R2)  A0=-A0 

SET  UP  THE  WINDOW  ALTITUDE  PLUS  ITS  UPPER  AND  LOWER  . 

3  CU  NEAR  I E  S . 

H  UP=  HNO  W 

HH=H NOW-WIN DOW/ 2. 

H  D  N  =  HNO  W  -  WT  N  EC  W 

SET  UP  CONTROL  INTEGERS  FOR  EACH  RAY 
DO  220  1=1, NELC 
IF  (H  (I,  K)  -HUP)  200,  200,210 
200  IF  (H  (I  ,  K)  -HEN)  21  5,205,205 
205  L (I) =2 

GO  TO  220 
2  10  L  (I)  =3 

GO  TO  220 
215  L (I) =1 
220  CONTINUE 

SET  UP  INITIAL  CONDITIONS 
E  =  0 . 0 


J  SLC?E=  0 
LC  =  C 
K  T  C  T  =  0 
:i  EL  =  NELO  -  1 
DO  225  I  F  S=  1 , 2 
TDIFF  (IFS)  =0.0 

225  ESIG  (IFS)  =0.0 
C 

C  SIT  0?  A  ICC  t  WHICH  SCANS  RAY  HEIGHT  VERSUS  EL  E  V  AT  I 

C  ANGIE. 

DO  6 CO  1  =  1,  Hr  L 
IF  (1-  1)  226,226,230 

226  IF  (1  (1)  - 2)  230,60  0,230 
230  J  LA  £T=J  SLC  EE 

JC=IASS  (I  (I)  -I  (1+1)  ) 

IF  (I  (I)  -I  (1  +  1)  )  260,600,240 
C 

C  NORMAL  RAY  ORDER 

240  J  SICPE= 1 

IF  (JLAST-  1)  280,280,31  0 

C 

C  INVERTED  RAY  ORDER 

260  J  SI CFE=  2 

IF  (JLAST-1)  280,3  1  0,280 
C 

C  EITHEP  ONE  OR  TWO  WINDOW  LIMIT  CROSSINGS 

280  IF  (JC-1)  30C,300,  380 

C 

C  ONE  WINDOW  LIMIT  CROSSING 

300  L C  =  LC  + 1 

IF  (LC-1)  34C,340,  360 
C 

C  RAY  ORDER  REVERSAL  OCCURS.  CHECK  IF  ONE  OR  TWO 

C  WINDOW  LIMIT  CROSSINGS. 


1  16 


310  IF  (JC-1)  320,3  20,230 


ONE  WINDOW  LIMIT  CROSSING.  CHECK  IF  RAY  REVERSAL 
OCCCRS  INSIDE  OR  0  NT  SIDE  THE  WINDOW 
320  IF  (LC)  2  30,230,560 

FIRST  WINDOW  LIMIT  CROSSING 
340  J1=I 

GO  TO  6  0  C 

SECOND  WINDOW  CROSSING 
36  0  J  2  =  1  +  1 

FINE  WHICH  RAYS  LIE  ABOVE  AND  BELOW  THE  WINDOW  CENTE 
DO  370  J12=J1,J2 
IF  (JSLOPE-1)  362,  362,  364 
362  IF  (HH-H  (G12+1,  K))  370,366,366 
364  IF  (HH-H  (J12  +  1,K))  366,366,370 

366  K 1 = J  12 
K  2=  J  12+  1 
GO  TC  440 
370  CONTINUE 


30  TH  WINDOW  LIMITS  CROSSED.  CHECK  FOR  "RADIO  HOLE" 
330  IF  (DABS  (H(I,K)  -K(I+1,K)  )-5.0*HT3)  390,390,  5o0 

BOTH  WINDOW  LIMITS  CROSSED.  NO  "RADIO  HOLE"  HAS  BE 
F  CU  ND. 

390  J 1=1 
J  2= 1  +  1 
K  1  =  J  1 
K  2  =  J  2 


EN 


THE  INC  R.  E  HEN  I  GE  RAYS  AT  THE  WINDOW  HAS  BEEN  ROUND. 
CHECK  IF  RAY  CEDZF.  IS  HOFZAL  OR  INTERRED. 

4  40  I?  (JSLCPE-1)  460,460,480 

NORMAL  RAY  CEDE  F.  CALCULATE  THE  UPPER  AND  LOWER 
ANGLES  SUETENDEZ  3Y  THE  KINDGF. 

46  0  E  UP  =  ELO  (J  1)  -EE N CT*  (H  ( J  1  ,  K)-HU?)  /  ( H  ( J  1 ,  K)  -  K  (J  1  +  1 ,  X ) ) 
EDN=ELO  ( J 2)  +EENCT*  (KDN-H  (J2 ,  K)  )  /  (H  ( J2- 1  ,  K)  -  H  ( J2,  K)  ) 
S  =  D  EQRT (EUP-EDN) 

?  SC  AI  =  P  H I  (J1  ,K) 

R SC AT=R HC  (J1  ,K) 

GO  TC  500 


INVERTED  RAY  CREER.  CALCULATE  THE  UPPER  AND  LOWER 
ANGLES  SUETENDEE  3Y  THE  WINDOW. 

480  E  UP  =  ELO  ( J  2 )  +DENCT*  (H  (J2,K)-HUP)  /  (H  (J2,  K)  -  H  ( J2-  1 ,  K) ) 
EDN  =  ELO  ( J  1 )  -EENCT*  (HDN-H  (J  1  ,K)  )  /  (H  ( J 1  +  1  ,  K)  -H  (J1,  K)  ) 
E  =  DSCF.r  (EEN-EOP) 

PSCAT  =  PH I  ( J 2  ,K) 

RSCAT=P. HC  ( J2  ,  K) 


SET  UP  INITIAL  CONDITIONS  TO  CALCULATE  MEAN 
ANGLE  AND  MEAN  TIME  CF  AFRIVAL. 

500  K  TGT=KTO  T  + 1 


ELEVATION 


TSIGT  (KTCT)=G(K1,K)  -  (G  (K1,K)-G(K2,K)  )  *(H(K1,K)-HH) 
V  (H  (K1,  K)  -H  ( K 2 ,  K)  ) 

A LP H=-S  (K1,K)+(S(K1,K)-S  (K2,  K)  )  *  (H  (K  1  ,  K)  -  HH) 

1/  (H  ( K 1 ,  K )  -H  (K2  ,  K) ) 


CALCULATE  ANGLE  INCREMENT  OF  EACH  WAVEFRONT 
ES  IGT  (KTCT)  =E*RSCAT*  C3QET  (DCOS  (  A0)  /DCCS  (ALPH)  ) 
PSCT  (KICI)  =  f  SCAT 


C  RESET  INITIAL  CONDITIONS  TO  CONTINUE  SCAN  OF 

C  HEIGHT  VS  ANCLE. 

56  0  LC  =  C 
600  CONTINUE 
C 

C  CHECK  IF  WINDOW  IS  IN  A  SHADOW  REGION  RESULTING  E  RON 

C  CAUSTIC. 

IF  (KTOT)  790,760,  6  10 
C 

C  WINDOW  IS  NCI  IN  A  S HA  DO R  REGION.  PUT  WAVEFRONTS  IN 

C  ORDER  OF  THEIR  TIM ES -OF- ARRIVAL . 

610  CONTINUE 

DO  660  I S=  1 , KTCT 
T  MI N  =  10  . 0 
DC  650  IR=1 ,KTCT 
T  MI  N  =DM  I N  1  ( T  M I N , TS IG  T  ( I R  )  ) 

IF  (TMIN-1SIGT  (IR)  )  640,630,  640 

630  I  M I N  =  I R 
640  CONTINUE 
650  CONTINUE 

TSIG  (15)  =  T  SI  GT  (-IMIN) 

ESIG  (IS)  =  E  S  I  GT  (IMIN) 

?  SC  (IS)  =  ESC  I  (IMIN) 

TSIGT  (IMIN)  =20.  C 
660  CONTINUE 
C 

C  CHECK  IF  MORE  THAN  ONE  WAVEFRONT  IS  PRESENT  IN  THE 

C  WINDOW. 

IF  (KTOT-  1)  700  ,700,665 

C 

C  CALCULATE  TIME  DIFFERENCE  OF  ARRIVAL  (X.D.O.A.) 

C  BETWEEN  WAVEFRONTS. 

665  CONTINUE 


1  1  9 


Ei 


c 

c 


c 

c 


c 

c 

c 


.10  6 71  11=2, KTCT 
c7  j  ::  IE  F  (IT)  =:  ill  (IT)  -1  3IG  (  1) 

5  21  UP  A  I C  C  £  1C  CALCULATE  TOTAL  FIELD  STRENGTH  Ch 
PCFEFE  TENSITY  I?  WA7EFECNTS  OVERLAP  IN  TINE. 

7  0  0  E  s  u  m  =  o .  J 

Du  721  I F 1=  1 ,  KTCT 
IF  (TDIFF  (IFI) -TPW)  710,710,  720 
7  11  ZSUM=E3UM+E3IG  (IFI)  *  LCDS  (FSC  (IF  I)  -2.  D*?I*  FZ.Q 
1*11  IFF  (IFI)  ) 

72  1  CONTINUE 

ZSUMOALS  (ESLM) 

CHECK  IF  FIELD  STRENGTH  IS  TOO  LON  FOR  CALCULATION 
IF  (E3UM/EDIFF-C.  0000000  1)  730,  730,730 

CHECK  IF  FIELD  STRENGTH  CF  PONES  DENSITY  IS  TO  3E 
CONFUTED. 

730  IF  (NPLCT-1)  740,750  ,760 
740  If  (KPLOI-  1)  750,750,760 

CALCULATE  TOTAL  RELATIVE  FIELD  STRENGTH. 

75  0  5?=10.*DICG10  (E SUM/ EDI F F) 

GO  IC  300 

CALCULATE  TOTAL  RELATIVE  POWER  DENSITY. 

7  c  0  3  P  =  20. *  DLCG 1 C  (E  SUM/EDI  ?  F ) 

GO  TC  300 

SET  FIELD  ST FE NOTH  AND  POWER  DENSITY  FOR  A  SHADOW 
R EG ICM. 

780  3P  =  -1 000CCO. 0 


120 


PRINT  HEADING  AND  RAY  HEIGHTS  IF  THIS  IS  A  NEW 
DISTANCE . 

800  IF  (J-1)  810,810,820 

810  CONTINUE 

XNAUT=VX  (K)/CNAUT 
IF  (KPLOT-1)  616,812,814 
812  PRINT  1200,XNAUT 
GO  TO  816 

814  PRINT  12  10,XNAUT 

POSITION  THE  PEN  AND  DRAW  THE  NEW  CR  DIM  AT  E  AXES  IF 
THERE  IS  A  PLOT. 

816  IF  (NPLOT )  820,820,8  18 
818  CALI  PLOT  (SNGL  (XN  )  ,  0.  0  ,  - 2) 

CALI  PLOT  (0.0,  SNG  L  {  YS)  ,3) 

CALI  PLOT  (0.0,  SNG  L  (  YN)  ,2) 

CALL  PLOT  (0.0, 0.0, 3) 

I  CO  DE=  3 


PRINT  THE  RELATIVE  FIELD  STRENGTH  OR  POWER  DENSITY  A 
NUMBER  OF  WAVEFRONTS  IN  HEIGHT-GAIN  CALCULATIONS. 

820  IF  (KPLOT)  840,640,830 
330  PRINT  1220,J,HH,S? 


SET  PARAMETEES  IF  THERE  IS  A  PLOT. 
840  IF  (NPLOT)  900,900,  850 
850  Y  =  HH/YDI V 


CHECK  IF  FIELD  STRENGTH  CF  POWER  DENSITY  IS  TOO  LOW 
FOR  PLOTTING. 

IF  (SP+40.0)  660,870,870 


c 


FIELE  STRENGTH  CE  POWER  TENSITY  LESS  In  AN  -4 9  D3. 


C  TURN  OP  F  ELCTTEr. 

8o0  X=-40.0/ZEIV 

CALI  PLOT  (SNGI(X)  ,SNGL(Y),ICOD£) 

I  CO  E  £=3 
GO  TC  900 
C 

C  PICT  TOTAL  RELATIVE  FI2LE  STRENGTH  CE  POWER  DE 

870  X  =  S  P/ZD  I V 

IF  (J-1)  £80,980,382 

880  I  CO  E  E=3 
GC  TC  884 
882  I  CO  E  £=2 

884  CALI  PICT  ( SNGL  (X)  ,  S  NG  L  {  Y)  ,  ICOD  E) 

9  00  CONTINUE 
1000  CONTINUE 
C 

C  POSITION  THE  PEN  IF  THERE  HAS  BEEN  A  PLOT. 

1100  IF  (NPLCT)  1  120,1120,1  1  10 
1110  CALL  PLOT  (0.,  0-,  —  999) 

CALL  PLCT  (2. 0,2.0,  -3) 

CALL  FACTOR  (SNGL  { SC  ALE)  } 

1 120  CONTINUE 
RETURN 
C 

1200  FORMAT  (///, IX, 'PRINTOUT  CF  FIELD  STRENGTH  AND 
I'T.D.O.A.  S  AT  X  =•  ,F8. 2, 2X, 'NAUTICAL  MILES', 
2//5X, 'NO. ',8X,' HEIGHT  ( FT) ’ , 6X  ,  ’ FI  ELD  STRENGTH 
1210  FORMAT  (///, IX, 'PRINTOUT  CF  PCWER  DENSITY  AND 
I’T.D.O.A.  S  AT  X  =•  ,F8. 2, 2X, 'NAUTICAL  MILES', 
2//5X, 'NO .' ,SX, ’ HEIGHT  ( F T) ', 7X ,' POWER  DENSITY 
1220  FORMAT  ( 4  X,  1 3 , 8  X,F  1  0 . 2  ,  1  OX,  E  1 2  .  2) 

END 


NSITY. 
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S  U3 ROUTINE  FLCTTR 
C 

C - this  routine  sets  up  the  axes  foe  ail  plots. 

c 

IMPLICIT  FLAI*8  (A-H,  O-Z) 

REA1+4  AHL,T1I,  YTL,YTL,T1E 

COMMON  /TWOCCM/  RA , DELX , AES, A 30 , ELC (5 1) , E XM AX, XFINAL 
1  ,  JR  A  Y 

COMMON  /THE COM/  VX  (30)  ,  YMIN,YMAX,XDIV,YDIV 

1  ,  H  (51,30)  ,0(51,30)  ,  H  N  (100)  ,  RN  (100)  ,  S  (51,30) 

2  ,JPLT,KEEF,KGRAD, KP LCT, JPL CT , IPLACE,  JCASE, NELO 

3  , NR E E, NRA Y, NPRO,N PLOT, NGRAD, SCALE 

COMMON  /FCRCCM/  DN  DZ  (1  00 )  ,R  M  (1 0  0)  ,  DMDZ  ( 10  0)  ,  KMI E 
1  ,  K  MGR AE,  NMIE  ,  NMGRAD,HF  (1  0  0)  ,  DATCTL 

DIMENSION  AHL  (18)  ,TTL(18)  ,XTL(1  8)  ,  YTL  (18)  ,TLE(18) 
DATA  YSIZ,  SIZE/5. 0D0,1. 505D-1/ 

DATA  IZTI , ZSIZ, ZMI N , ZDI V/  1 0, 1 . 0  DO 
1  ,  -  2  .  OD 1 ,  4  . CD  1/ 

C  SET  UP  INITIAL  CONDITIONS  IF  THIS  IS  THE  FIRST  PLOT 

IF  (JPL  T )  13,10,20 

10  YDI  V=  (YMAX-YKIN)/YSIZ 
YDI VA=YDI V/1  COO.O 
Y MI NA=Y MIN/ 1000.0 
YPT  =  YSIZ*1.3 
ZPT  =  YSIZ  +  0.  5 
YS= (AHS- YMINJ/YDI7 
J?LT=1 
C 
C 


READ  PLOT  SCALE  FACTOR  ANI  PLOT  BANNER. 
REAL  110, SCALE 

REAL  120 ,  IAHL,  (AMD (K) , K=  1,18) 


c 


C  INITIALIZE  TEE  HOT  TEE  A  ND  SCALE  THE  PLOT  SIZE. 

CALI  PLOTS  (0,0,0) 

CALL  PLOT  (0.,0  .,-999) 

CALI  SET  M  SG  (2) 

CALL  FAC  TOE  (SNGL  (SCALE)  ) 

C 

C  POSITION  TEE  PEN  AND  PRINT  THE  FLCT  BANNER. 

15  CALL  PLOT  (S. 0,8. 0,-3) 

CALL  SYMBOL  (0. 0,0.0  ,SN3I  (SIZE)  ,AHD,00.  0,  IAHD) 

C 

C  POSITION  THE  PEN  FOR  THE  NEXT  PLOT. 

CALI  PLOT  (0.  ,  0  . ,-  999) 

CALI  PLOT  (2. 0,2. 0,-3) 

CALL  FACTOR  (SNGL  (SC  ALE)  ) 

C 

C  READ  AXIS  IAEE1S  AND  TWO  IINES  CF  PLOT  TITLES. 

20  READ  120  ,  IXT1,  i  XTL  ( K)  ,  K=  1 ,  1  8) 

READ  120,IYTI,  (YTL(K)  ,K=1,18) 

READ  120,  ITTL,  (TTL  (K)  ,  K=  1,  18) 

READ  120,  ITLE,  (TLE  (K)  ,  K=  1,18) 

C 

C  CHECK  WHICH  ELCT  THIS  IS. 

IF  ( JPLQT-5)  22,24,26 
22  CONTINUE 

IF  (JPLOT-2)  30,40,60 
C 

C  SET  UP  INITIAL  CONDITIONS  FOR  MODIFIED  INDEX  OF 

C  REFRACTIVITY  PROFILE  AND  GRADIENT  PLOTS. 

24  CONTINUE 

XMAX  =  520.  0 
XMIN  =  0.0 
GO  TO  50 
26  CONTINUE 

124 


XMAX  =  350.0 
X  MI N  =  0.0 
GO  TC  50 
C 

C  SET  OP  INITIAL  CONDITIONS  FOR  R  EFE  ACT  17  IT  Y  PF.OFIL 

C  AND  GRADIENT  PLCTS. 

30  X  MA  X=400 . 0 
X  MI N =0. 0 
GO  TC  50 
40  XMAX  =20  0.0 
X  MI N=-6  0  0.0 
50  X SI Z=4. 0 

XDIV=  (XMAX-XMIN)  /X  SI  Z 
XDIVA=XDIV 
GO  TC  70 
C 

C  SET  OP  INITIAL  CONDITIONS  FOR  RAY  TRACE,  FIELD 

C  STRENGTH,  AND  PCWER  DENSITY  PLOTS. 

60  XNAX=EXMAX 
X  M I N  =  0  .  0 
XSIZ=10. 0 

X  DIV= (XMAX-XMIN)/XSIZ 
XDIVA=XFINAL/XSIZ 
XEE=XSIZ-C. 5 
C 

C  SET  UP  PARAMETERS  FOR  AXIS  LABELS  AND  PLOT  TITLES 

70  XTTL=  (XSIZ-SIZE*I7TL)/2. 

XTLE=  (XSIZ-SIZE*ITLE)/2. 

C 

C  DRAW  AND  IABEI  THE  AXES. 

CALL  AXIS  (0.0,  0.0,  YTL,  IYTL,  SNG  L  (YSIZ)  ,90.0, 

1SNGL  (YMINA)  ,SNGI(YDIVA)  ) 

CALI  AXIS  (0.0,0.0,XTL,-IXTL,SNGL(XSIZ)  ,0.0, 
ISNGL(XMIN)  ,  SNG  I  (XDIVA)  ) 


125 


c 


C  PRINT  THE  FLC1  TITLES. 

74  CALL  SYMECL  (S  NGL  { XTTL)  ,  SNGL  (Y?  T)  ,  SNGL  (SI  ZE)  ,  T7L,  0.  0  , 
II ITL) 

CALL  SYMECL  (S NGL  (XTLE)  ,  SNGL  (ZP  T)  ,  SNGL  (  SI ZS)  ,  TLE,  0 . 0  , 
1ITLE) 

C 

C  DRAW  THE  REFERENCE  SCALE  FOR  ALL  FIELD  STRENGTH  AND  . 

C  POWER  DENSITY  PLOTS 

IF  (JPLOT-5)  78,82,  82 
78  CONTINUE 

IF  (JPLGT-3)  82,32,80 

80  CALL  AXIS  (SNGL  (XRE)  ,SNG1  (ZPT)  ,  10H  GAIN  (D3),-IZTL, 

IS  NGL  (ZSIZ)  ,0.0,SNGL(ZMIN)  ,SNGL(ZDIV)  ) 

CALI  PLOT  (SNGL  (XSIZ)  ,SNGI  (ZPT)  ,3) 

CALI  PLOT  (SNGL  (XSIZ)  ,SNGI  (YPT)  ,2) 

82  CALL  PLOT  (0.0, 0.0, 3) 

C 

C  DRAB  THE  EARTH  SURFACE  IF  NOT  AT  ZERO  MEAN  SEA  LEVEL. 

IF  (AHS-YMIN)  86,36,84 
84  CALL  PLOT  (0 . 0  ,  SNG  L  (  YS)  ,  3) 

CALL  PLOT  (SNGL  (XSIZ)  ,SNGI  (YS)  ,  2) 

CALL  PLOT  (0.0, 0.0, 3) 

C 

C  POSITION  THE  PEN  FOR  THE  REFR ACTIVITY  GRADIENT 

C  PROFILE  ELCT. 

86  CONTINUE 

IF  (JPLCT-2)  88,90,88 
38  CONTINUE 

IF  (JPLCT-5)  1  00,90,90 
90  X 11=  -XaiB/XDIV 

CALI  PLOT  (SNGL  (XN)  ,0.0, -3) 


100  CONTINUE 
105  RETURN 


c 


1  10  FORMAT  (3F15.7) 

120  FORM  AT  ( I  2 , 8  X,  1  7A4  ,  A  2) 
£  ND 


C 


Q  ****************************************************  * 


BLOCK  DATA 

C  SUBROUTINE  BLOCK  DATA. 

IMPLICIT  REAL*8  (A-H,  O-Z) 

COMMON  /CNECCM/  CRH , CE X , CFG , DTE , PK AS E , T H£TA , C, C 1 , C2 

1  ,  R  HC  (51,30)  /  PHI  {51,30}  ,  ?I,  CMF,CNAUT,  DENOT 

2  , FR  Q ,TE  W, ABSRH, NHV,  NSL, NEM S, JARS A, NDATA 
COMMON  /THOCCM/  RA  ,  DELX  ,  AHS  ,  AHO  ,  ELC  (5  1)  ,  S  XM  AX  ,  XFI NAL 

1  ,JRAY 

CCMMCN  /THRCCM/  VX  (3  0)  ,  Y  MIN,  YI1A  X  ,  XDI V  ,  YDI  V 

1  ,H  (51,30)  ,G  (5  1,30)  ,  EN  (100)  ,  RN  (100)  ,3  (31,30) 

2  ,JPLT,KFEF,KGRAD, KP ICT, JPLOT, IPLACE, JCASZ, NELO 

3  , NREF,  NEAY,NPRO,NPLCT,NGE AD, SCALE 
DATA  R HO,  EH  1/1  530*  0.  DO  ,  1  530*0.  DO/ 

DATA  ELO/  5  1  *  G.  DO/ 

DATA  VX,  E  ,G  ,  H  N  ,  FN/3  0*0.  DO, 15  30  *0. DO,  1  53  0*  0  .  DO  , 
1100*0. DO,  100*0. CO/ 

DATA  3/  1  5  30  * C.  D 0/ 

DATA  DTE, FA, FI/ 1.7 45  3D-0  2,2. 092  5D+07,3.  14  15  9D0/ 

DATA  CMF,CNALT/3.28  1  DO,  6.G76D3/ 

END 


/* 

//GC.SYSIN  DD  * 

1 

0 

2 

8 

16250.  36.330 
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2672.  233. 4 C 1 

2595.  236. C51 

1480.  274.466 

14o6 .  272.235 

333.  3  1  4. S85 

72.  35  1.  813 

23.  357.481 

75.  45  90. 

70.  00  5 0  0  C. 

20.  00  200.00 

0.50  -0.50 

2900.  00  6.50 

1  2  3 

1  1  C  0  1  1 

111111 

1.00 

12  19  APB  1584 

22  REFRACTIVITY  (N-UNITS) 

38  HEIGHT  AEOVE  EARTH  {THOUSANDS  OF  FEET) 

32  ATMOSEHEBIC  PEFR  ACTIVITY  PROFILE 

42  EHAHRAN,  SAUDI  ARABIA  -  1  MAY  1978,  0000Z 

34  REFRACTIVITY  GRADIENT  (N-UNITS/KM) 

38  HEIGHT  AEOVE  EARTH  (THOUSANDS  OF  FEET) 

29  REFRACTIVITY  GRADIENT  PROFILE 

42  EHAHRAN,  SAUDI  ARABIA  -  1  MAY  1978,  0000Z 

40  MODIFIED  INDEX  OF  REFRACTIVITY  (M-UNITS) 

38  HEIGHT  AEOVE  EARTH  {THOUSANDS  OF  FEET) 

38  MODIFIED  INDEX  OF  REFRACTIVITY  PROFILE 

42  EHAHRAN,  SAUDI  ARABIA  -  1  MAY  1978,  OOOOZ 

42  MOD  INDEX  OF  REFRAC  GRADIENT  (M-UNITS/KM) 

38  HEIGHT  AEOVE  EARTH  (THOUSANDS  CF  FEET) 

43  MOD  INDEX  OF  REFRACTIVITY  GRADIENT  PROFILE 

42  EHAHRAN,  SAUDI  ARABIA  -  1  MAY  1978,  OOOOZ 
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37  El  STANCE  ALCN3  EARTH  (HA  DTI  CAL  MILES) 

38  HEIGHT  AECVE  EARTH  (THOUSANDS  OF  FEET) 

56  FAY  TEACF  FOR  DHAHRAN,  SAUDI  AF.AEIA  - 

♦LINE  CCN'T*  1  MAY  73,  0000Z 
56  HO  =  SO  FT  RAY  ANGLES  FROM  -0.5  TO  0.50  .. 

♦  LINE  CCN»  T*  DEGRESS 

37  DISTANCE  ALONG  EARTH  (NAUTICAL  MILES) 

38  HEIGHT  AECVE  EARTH  (THOUSANDS  OF  FEET) 

58  EEL  ELD  SIR  -  DH  AHR  A  N, S  A UDI  A  ARABIA  - 

♦LINE  CCN'T*  1  MAY  78,  0000Z 

56  HO  =  SO  FT  RAY  ANGLES  FROM  -0.5  TO  0.50  .. 

♦LINE  CCN  *  T*  D  EG  E  EES 

Q  ************************************  ****************  ******* 
Q* *********************************************************  * 

EXAMPLE  OF  THE  OTHER  LATA  SET: 

1 

1 

2 

9 


16550.  000 

38.939 

2637  .2  87 

-6.  760 

150.  238 

3161. 000 

20  3.  4S3 

699.770 

-19.  978 

137. 021 

2478. 00C 

217.138 

606.184 

-5. 109 

151.  892 

2209. 000 

218. 512 

565.325 

-30.  442 

126. 558 

1488. 00C 

240. 461 

474.077 

-41.  159 

115. 841 

627. 000 

275. 8S9 

374.338 

-34. 524 

122. 477 

374. 000 

284.634 

343.35 1 

-164.  929 

-7.  929 

49.  000 

338. 235 

345.928 

26.  602 

183.603 

23. 000 

337.544 

341.155 

200. 000 

350. 000 

75.  45 

SO. 

70. 00 

5000. 

20.  00 

2CC.0C 

0.  50 

-0.  50 

9800.00  1340.00 


1 

1 


12 

22 

38 

32 

42 

34 

38 

29 


4  2 


40 

38 

38 

42 

42 
38 

43 
42 

37 

38 
56 


56 


37 

38 
58 

56 


C****** 


C  C  0  0  0 

112  11 

1.00 

26  AFF  198  4 

REF P.  AC  II VII  Y  (N -  UNITS) 

HEIGHT  AEOVE  EAF.IH  (THOUSANDS  OF  FEET) 
ATMOSPHERIC  F.E F F.  AC T I  V IT Y  PROFILE 
EHAHRAN,  SAUDI  ARABIA  -  23  MAY  1978,  OOOOZ 
F.EFR  ACTIVITY  GRADIENT  (N-UNITS/KM) 

HEIGHT  AEOVE  EARTH  (THOUSANDS  CF  FEET) 
REFRACUVITY  GRADIENT  PROFILE 
DHAHFAN,  SAUDI  ARABIA  -  23  MAY  1973,  OOOOZ 
MODIFIED  INDEX  OF  REFRACUVITY  (M-UNITS) 

HEIGHT  AEOVE  EARTH  (THOUSANDS  OF  FEET) 

MODIFIED  INDEX  OF  REFRACUVITY  PROFILE 
EHAHRAN,  SAUDI  ARABIA  -  23  MAY  1978,  00002 
MOD  INDEX  OF  REFRAC  GRADIENT  (M-UNITS/KM) 
HEIGHT  AEOVE  EARTH  (THOUSANDS  OF  FEET) 

MOD  INDEX  OF  REFRACUVITY  GRADIENT  PROFILE 
DHAHRAN,  SAUDI  ARABIA  -  23  MAY  1978,  0000Z 
DISTANCE  AICNG  EARTH  (NAUTICAL  MILES) 

HEIGHT  AEOVE  EARTH  (THOUSANDS  OF  FEET) 

RAY  TRACE  FOR  DHAHRAN,  SAUDI  ARABIA  - 

♦  LINE  CCN’  T*  23  MAY  78,  OOOOZ 

HO  =  SO  FT  RAY  ANGLES  FROM  -0.5  TO  0.50 

♦LINE  CCN  *  T*  DEGREES 

DISTANCE  ALONG  EARTH  (NAUTICAL  MILES) 

HEIGHT  AEOVE  EARTH  (THOUSANDS  OF  FEET) 

SEL  FID  STR  -  DH AHR A N,S AUDI  A  ARABIA  - 
♦LINE  CCN 1 T*  23  MAY  78,  OOOOZ 
HO  =  90  FT  RAY  ANGLES  FROM  -0.5  TO  0.50 

♦  LINE  CCN1!*  DEGREES 
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